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THEORY OF SEPARATION OF VARIABLES FOR LINEAR PARTIAL DIFFERENTIAL 
EQUATIONS OF THE SECOND ORDER IN TWO INDEPENDENT VARIABLES 

by Marvin E. Goldstein 
Lewis Research Center 

SUMMARY 

Necessary and sufficient conditions which any linear second-order partial differential 
equation in two independent variables must meet whenever it can be transformed into a 
separable equation are given. These conditions are used to develop a calculational' pro- 
cedure for determining whether any given equation of this type can be transformed into a 
separable equation and also to develop a procedure for determining the various changes 
of variable which will lead to separable equations. 


INTRODUCTION 

Perhaps the most useful way of obtaining solutions to linear partial differential equa- 
tions is the method of separation of variables. Unfortunately this method is only appli- 
cable to a small number of equations. However, the applicability of this method can be 
increased somewhat if the variables in the differential equation are transformed before 
the method is applied. It is therefore useful to be able to tell whether a given partial dif- 
ferential equation can be transformed into a separable equation (that is, an equation which 
can be solved by the method of separation of variables) by changing both its dependent 
and independent variables. 

We shall therefore give necessary and sufficient conditions which the coefficients of 
any second-order linear partial differential equation in two independent variables must 
satisfy in order that the equation be transformable into a separable equation. These con- 
ditions require that the coefficients (or more precisely, certain combinations of the coef- 
ficients) be expressible in certain functional forms. Since it may not always be easy to 
tell in practice simply by inspection whether a given set of coefficients can be expressed 
in this way, alternate forms of these conditions are given which allow the coefficients of 
a given equation to be tested by direct calculation. In order to use the procedures de- 



veloped for this purpose, it is required in the worst situation that the solution to an or- 
dinary differential equation be found. 

In addition, if it is found that a given equation can be transformed into a separable 
equation, then the formulas developed herein can be used to calculate all the possible 
transformations which will bring the equation into separable form. 

Some incomplete or limited studies along these lines have already been carried out. 
Thus, in references 1 and 2, all the possible conformal transformations which trans- 
form the separable equation 


V 2 U + k 2 U = 0 

o 

into another separable equation have been enumerated. (V will always be used herein 

2 

to denote the two-dimensional laplacian. ) In reference 1, k is taken as a constant and, 
2 2 

in reference 2, k is taken to be a constant divided by y where y is one of the orig- 
inal independent variables. In reference 3, sufficient conditions for the equation 

= xu 

ax ay 

to be transformable into a certain type of separable equation are given. Boussinesq 
(ref. 4) showed that the equation 


V 2 U +Vcp • VU = 0 

where <p is a given harmonic function, can always be transformed into a separable equa- 
tion by transforming both the dependent and independent variables, the tranformation of 
the independent variables being a conformal transformation. A slight generalization of 
Boussinesq's result is given in reference 5. It was shown there that the potential could 
be a slightly more general function. 

All these results will emerge as special cases of the general theory developed here- 
in. Limited results for special equations in more than two independent variables are 
given in references 1 and 5 to 9 . 

We begin by finding the restrictions that are imposed on an equation by the require- 
ment that it be separable. It turns out that it is convenient to distinguish between those 
equations for which the method of separation of variables leads to two ordinary differential 
equations of the highest possible order consistent with the type of partial differential equa- 
tion and those for which it does not. Only the former type of equations are called separ- 
able herein. The latter type are called ’’weakly separable but not separable.” However, 
since the method of separation of variables often leads to useful results even when an 


2 


II 


equation is weakly separable but not separable (ref. 10), this case is also considered in 
detail. 

The various algebraic transformations which can be applied to a differential equa- 
tion are discussed. It is then determined what restrictions are imposed on these trans- 
formations by the requirement that they transform a given partial differential equation 
into a separable equation. Once this is done, the restrictions imposed on the coefficients 
of the original equation can be found. It turns out that this is best accomplished by con- 
sidering the elliptic, parabolic, and hyperbolic equations separately. 

We shall assume that all the functions which are encountered can be differentiated 
as many times as is necessary. We shall say that the function f of two variables is not 
equal to zero and write f ^ 0 if it takes on the value zero only at discrete points (or, at 
most, along line segments). 


SEPARABLE EQUATIONS 


The most general second-order linear and homogeneous partial differential equation 
in two independent variables has the form 


aifH.2,3 

9£ 2 


3 2 U 

3 | dr] 




ar 


3£ 


31 


( 1 ) 


where the coefficients a, 0, y, A, B, and C are real functions of | and ??. We shall 
suppose that a, 0, and y are not all zero, for this would imply that equation (1) was in 
reality a first-order equation. We shall also require that a, 0, and A (or y, 0, and 
B) are not all zeros. Otherwise, equation (1) could be essentially an ordinary differen- 
tial equation. 

2 

Equations of this type are further classified by the sign of the discriminant 0 - ay. 
Thus, if 


0 2 - ay > 0 (2) 

the equation is said to by hyperbolic . If 

0 2 - ay = 0 (3) 

the equation is said to be parabolic , and if 


3 



flP - ay < 0 


(4) 


the equation is said to be elliptic . We shall always assume that the domain of definition 
of the equation has been restricted in such a way that the sign of the discriminant does 
not change. This assumption will simplify the following presentation but will not affect 
its generality. It is well known that the type of boundary value problems which can be 
solved by an equation of the form (1) depends only upon whether this equation is hyper- 
bolic, parabolic or elliptic. 

When the coefficients of equation (1) satisfy certain restrictions, this equation can 
be solved by the method of separation of variables. This method consists of substituting 
the trial solution 


m,n) = E(£)H(7?) 


(5) 


into equation (1) to obtain 

— (aS" + AS') +— (yH M +BHM +C +— 2/3E r H' = 0 (6) 

EH EH 

where the primes denote differentiation of the functions with respect to their arguments. 
Now suppose it is possible after division by (H , /H) n (H’/S) m f for some nonzero function 
f of £ and rj and n, m = 0 or 1, to write equation (6) as the sum of two terms, one of 
which is a function of £ only and the other a function of rj only, for all choices of the 
functions E and H. Then since £ and t] are independent variables, we can conclude 
that the trial solution (5) satisfies equation (1) only if each of these terms is equal to a 
constant (called the separation constant). This, in turn, implies (in view of the assump- 
tions made about the vanishing of the coefficients) that H and H must each satisfy an 
ordinary differential equation and that, if E and H do satisfy these equations, then 
equation (5) is indeed a solution of equation (1). Each of these ordinary differential equa- 
tions is of, at most, second order. Hence, S and H can each involve two arbitrary 
constants of integration. Since the separation constant is also arbitrary, it is clear 
that the solution (5) can contain at most four arbitrary constants. Thus, if the separ- 
ation of variables method works, it will lead to an, at most, four-fold infinite family of 
solutions to equation (1). If the family of solutions obtained by this method is suffi- 
ciently large, it is possible to express any reasonable solutions as a linear combination 
of members of this family. The family is then said to be complete. If equation (1) is 
either hyperbolic or elliptic, it is possible in some cases to obtain two second-order 
ordinary differential equations when the method is applied. 

Equation (1) will be called separable only if the method of separation of variables can 
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lead to two ordinary differential equations of the highest possible order consistent with 
the type of partial differential equation. Thus, if equation (1) is either hyperbolic or 
elliptic, it is said to be separable only if the method of separation of variables leads to 
two second-order ordinary differential equations. It will be shown below that in the par- 
abolic case the method can lead to, at most , one first-order and one second-order ordin- 
ary differential equation. Therefore, a parabolic equation will be called separable if the 
method of separation of variables leads to one first-order and one second-order ordinary 
differential equation. 

The method of separation of variables is sometimes useful (see ref. 10) even when 
it does not lead to ordinary differential equations of the highest possible order. There- 
fore, we shall call equation (1) weakly separable if the method leads to two ordinary dif- 
ferential equations regardless of their order. Notice that separation solutions to weakly 
separable equations involve at least two arbitrary constants and that every separable 
equation is weakly separable. 

In order that equation (1) be weakly separable, it is first necessary that, for arbi- 
trary S and H, equation (6) can be written as the sum of two terms one of which is a 
function of £ only and the other a function of r] only. It is clear that this cannot occur 
unless at least one of the coefficients a, 0, or y is equal to zero. 

First, suppose that equation (1) is elliptic, then condition (4) implies that a * 0 and 
y ± 0. This implies that equation (1) is both elliptic and weakly separable only if 0 = 0. 

In addition, the coefficients of S TI and H ,T in equation (6) never vanish and therefore 
the method of separation of variables, if it works, will always lead to two second-order 
ordinary differential equations! Thus * an elliptic aquation is Separable if, and only if, 
it is weakly separable. 

Next, guppsjSe that equation (1) is parabolic, Condition (3) shows that, if ]3 ^ 0, then 
Ot # 0 and y £ Os Hence, when equation (1) ii parabolic and weakly separable, we can 
conclude that ( since One of these coefficients must be sgro) 0 - 0, But H 0 = 0, condi- 
tion (3) shows that either a = 0 or > - 0, Thus* one of the second derivatives must be 
missing from equation (8), This pixjves that in the parabolic Case the method Of Separa^ 
tion of variables can lead* at most ^ to one first-Order and one SeCpnd=Prder ordinary 
differential equation, j£ equation (i) were weakly Separable but not separable! the pre= 
ceding remarks show that We would have a m 0 = y = 0, But this is Contrary to hypoth^ 
esis, HenCe! we conclude that a parabolic equation is Separable if 3 and only M 3 It Is 
weakly separable. 

Finally , suppose that equation (1) is hyperbolic , tn this case , equation (1) can still 
be weakly Separable evgn if 0^0, However, if 0 # 0 t then at least one of the COeffi= 
Cients a and y must be Zero if equation (1) is to be weakly separable, Thus^ at least 
one of the Second derivatives will not occur In equation (8)| and, therefore, the method 
of Separation Of variables cannot lead to two iecond-order Ordinary differential equations. 


§ 



This shows that equation (1) cannot be hyperbolic and separable unless [3 = 0. If (3 were 
zero, then condition (2) shows that a * 0 and y # 0. Hence, if equation (1) is hyperbolic 
and (3 = 0, then it is weakly separable if, and only if, it is separable. 

The preceding discussion allows us to arrive at the following conclusions: 

(Cl) In all cases equation (1) is separable only if 0 = 0 . 

(C2) Equation (1) is weakly separable but not separable only if it is hyperbo li c and 
0 * 0 - 

(C3) If equation (1) is parabolic and separable or if it is weakly separable but not 
separable, then either a = 0 or y = 0 . 

Now suppose that /3 = 0. Then equation (6) becomes 

— (aH" +AS’) +— (yH" +BH’) + C = 0 (7) 

S H 

It is clear that, for arbitrary S and H, this equation can be written as the sum of two 
terms , one of which depends on £ and the other only on rj if, and only if , the coeffi- 
cients a, y, A, B, and C can be expressed in the following forms: 


a(l,n) =f(5»»7) <4(0 

(8) 

rU,v) =f (!,»?) e 1 ( 7 ?) 

(9) 

A(l,?7) = f(S,T7) <4(0 

(10) 

B(|,7?) = f(Z,v) e 2 (? 7 ) 

(11) 

=f(S,77)[d 3 (0 + e 3 ( 7 l)J 

(12) 


where f * 0, d^ and dg are not both zero, e^ and eg are not both zero, and dj and 
e ^ are not both zero. These restrictions follow from the restrictions placed on the van- 
ishing of the coefficients of equation (1). Thus, the conditions (8) to (12), together with 
the condition 


0 = 0 (13) 

imply that equation (7) is weakly separable, and conclusion (C2) shows that they also impl. 
that equation (1) is separable. 

Conversely, suppose that equation (1) is separable. Then conclusion (Cl) shows that 
condition (13) holds. Thus, equation (7) must be expressible as the sum of two terms one 
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of which depends only on £ and the other only 77. But this implies that conditions (8) 
to (12) hold. Hence, we arrive at the following conclusion: 

(C4) Equation (1) is separable if, and only if, its coefficients satisfy conditions (8) 
to (13) . 

There are additional restrictions imposed on the functions d^ and e^ by the sign of the 
discriminant /3 2 - ay of equation (1). The conditions (8), (9), and (13) show that if equa- 
tion (1) is separable, then 


/3 2 - ay = -f 2 d-,e 


1 C 1 


(14) 


Suppose first that equation (1) is hyperbolic. Then equations (2) and (14) show that 

d 1 (^)e 1 (t?) < 0 

at each point (£,77) of the domain of equation (1) This shows that the sign of dj(|) is 
different from the sign of e^(r]) at each point. Since £ and 77 are independent varia- 
bles, this, in turn, implies either that 


or that 


d^ > 0 and e^ < 0 


dj < 0 and e-^ > 0 

However, in view of the symmetry of equation (1) and of conditions (8) to (13), no gener- 
ality will be lost if we assume that the first of these always holds . 

Next suppose that equation (1) is parabolic. Then equations (3) and (14) show that 


d 1 U)e 1 (77) = 0 

at each point (£,77) of the domain of equation (1) . Again since £ and 77 are independent 
variables, this shows that either 


e l 


= 0 


or 


d 


1 


= 0 
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We have already indicated that both these conditions cannot hold simultaneously. Hence, 
in view of the symmetry of equation (1) and of conditions ( 8 ) to (13), no generality will be 
lost if we assume 


e = 0 and d^ * 0 

Finally, suppose that equation (1) is elliptic. Then equations (4) and (14) show that 

d^Oejfo) > 0 

at each point (i;,? 7 ) of the domain of equation (1) This shows that the sign of d^(l) is 
the same as the sign of e^(r]) at each point (1, 77 ). Since £ and 77 are independent var- 
iables, this, in turn, implies either that 

dj > 0 and e^ > 0 


or that 


dj < 0 and e 1 < 0 

It can be seen, however, from conditions ( 8 ) to (13) that no generality will be lost if we 
assume that a minus sign has been absorbed into the function f. We therefore assume 
that the first of these conditions always holds. 

We have therefore shown that the functions d^ and e^ satisfy the following restric- 
tions : 

dj > 0 and e^ < 0 (15) 

if equation ( 1 ) is hyperbolic , 

dj # 0 and e-^ = 0 (16) 

if equation ( 1 ) is parabolic , and 

d x >0 and e x > 0 (17) 


if equation ( 1 ) is elliptic . 

Now suppose that 0 * 0, y * 0, and a = 0. Then equation ( 6 ) becomes 


8 



( 18 ) 


— A — +-L (yH” + BIT +CH) +— 2/3E* = 0 
E H' H’ E 


Since y and /3 are not zero, it is clear that, for arbitrary E and H, this equation 
can be written as the sum of two terms, one of which depends only on £ and the other 
only on 77 if and only if the coefficients a, /3, y, A, B, and C can be expressed in the 
following form: 


0 

II 

1 ? 

(19a) 

/3(I,t?) =f(?,77)d 1 (|) 

(20a) 

y(b,n) =f(l,T7)e 1 (77) 

(21a) 

A(£,t 7) = 0 

(22a) 

B(£,t?) = f(£,77)[e 2 (77) +d 2 (£)J 

(23a) 

3 

^co 

II 

1 ? 

O 

(24a) 


where f ^ 0 and d^ + 0. Similarly, if j3 * 0, a * 0, and y = 0, the resulting form of 
equation (6) can, for arbitrary S and H, be written as the sum of two terms with one 
of them depending only on £ and the other only on 77 if and only if the coefficients a , 
/3, y, A, B, and C can be expressed in the form 


«(!, 77) = f(^,7?)d 1 (^) 

(19b) 

/3(|,77) = f(^,77)e 1 (77) 

(20b) 

y(£,*?) = 0 

(21b) 

A(!,T?) = f(^, 77 )|d 2 (|) +e 2 ( 77 )] 

(22b) 

B(£,t 7 ) = 0 

(23b) 

C(I,t?) =f(|,77)d 3 (^) 

(24b) 


where f ^ 0 and e^ # 0. Finally, if /3 # 0, a = 0, and y = 0, the resulting form of 
equation (6) can, for arbitrary E and H, be written as the sum of two terms with one 
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of them depending only on £ and the other only on tj if, and only if, the coefficients a, 
P, y, A, B, and C can be expressed either in the form (19a) to (24a) or in the form 
(19b) to (24b). Thus, either the conditions (19a) to (24a) or the conditions (19b) to (24b) 
taken together with the condition 


p * 0 (25) 

imply that equation (1) is weakly separable and conclusion (Cl) shows that they also imply 
that equation (1) is not separable. 

Conversely, suppose that equation (1) is weakly separable but not separable. Then 
conclusion (C2) shows that condition (25) holds, and conclusion (C3) shows that either 
a = 0 or y - 0. Thus, the appropriate form of equation (6) (depending on whether a = 0 
or p = 0) must be expressible as the sum of two terms with one of them depending on ij 
only and the other on rj only. But this implies that either conditions (19a) to (24a) or 
conditions (19b) to (24b) hold. This shows that the following conclusion holds: 

(C5) Equation (1) is weakly s e parable but not separable if, an d only if, its coeffi- 
cients satisfy condition (25) an d either conditions (1 9a) to (24a) or conditions (19b) t o (24b) . 

The conditions obtained above show that only a very small percentage of all the 
second-order linear partial differential equations in two independent variables are even 
weakly separable. However, a somewhat larger percentage of the second-order linear 
partial differential equations can be transformed into weakly separable equations by 
changing either their dependent or independent variables. Hence, the usefulness of the 
method of separation of variables can be extended by using it in conjunction with a change 
of variables. We therefore develop a procedure for determining whether a given second- 
order linear partial differential equation can be transformed (by changing either the de- 
pendent or independent variables) into a weakly separable equation. In addition, when- 
ever a given equation can be so transformed, a method for calculating the appropriate 
change of variables will be given. 


ALGEBRAIC TRANSFORMATIONS OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS 

Only a change in the dependent variable which is of the form 

V(|,rj) =X(|, v)m,V) (26) 

where A is any nonzero function of £ 77 , will transform the linear homogeneous differ- 
ential equation (1) into another linear homogeneous equation. Hence, only transforma- 
tions of this type are appropriate for our purposes. The allowable transformations of the 
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independent variables are much less restricted. Any change in the independent variables 
of the form 


x =x(!,? 7 ) 

y = yU,Tj) 


where x and y are any functions of | and 77 such that 


d(x,y) + o 

9 ( 1 , 77 ) 


(27) 


(28) 


will transform equation ( 1 ) into another linear homogeneous partial differential equation. 
Hence, all transformations of the independent variable which satisfy condition (28) will 
be considered. 

It is important to notice (refs. 3 and 11) that both a change of variable of the type (26) 
and one of the type (27) will leave the sign of the discriminant invariant. Thus, any 
change of variable which is of interest for the present purpose will always transform 
hyperbolic equations into hyperbolic equations, elliptic equations into elliptic equations, 
etc. This fact is very useful in the following analysis. 

Up to this point it has been convenient to consider any two differential equations which 
have different coefficients as being completely different equations. We shall sometimes 
change this point of view slightly and consider two differential equations which can be 
transformed into one another by a transformation of the type (26) or of the type (27) as 
being different forms of the same equation. 

Recall now that every second-order linear and homogeneous partial differential equa- 
tion can, by a change of variable of the type (27) be transformed into one and only one of 
three canonical forms, depending on the sign of the discriminant. Thus, every 
hyperbolic equation can be put in the form (ref. 11) 



a 2 u au 

* * 



+ cU = 0 


(29) 


Every parabolic equation can be put in the form 


ifu+alH+biM 

3x 2 a* ay 


+ cU = 0 


(30) 


and every elliptic equation can be put in the form 
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(31) 


^u + i!u +a au +b au 

ax 2 a y 2 0X 0 y 


+ cU = 0 


where the coefficients a, b, and c can be any functions of x and y. 

Hence, no generality will be lost if we assume that this transformation has already 
been carried out and therefore that the differential equation under investigation is already 
in one of these three forms. We also see from the preceding remarks that even though 
it is necessary to specify four functions (since eq. (1) can always be divided through by 
one of the coefficients of its highest derivatives which cannot all be zero) in order to 
characterize any second-order linear homogeneous partial differential equation in two 
independent variables, at most only three functions (namely, the three coefficients a, b, 
and c appearing in the canonical form of the equation) need be known to determine 
whether the equation can be transformed into a weakly separable or a separable equation. 
We shall see that in fact only two functions need be known for this purpose. 

If the coefficient b in equation (30) is zero, this equation is essentially an ordinary 
differential equation and can be solved as such. We therefore exclude this case by im- 
posing the restriction b^O. With this restriction it is easy, though somewhat tedious, 
to verify that equations (29) to (31) can never be transformed by a change of variable of 
the type (26) or of the type (27) into an equation whose coefficients violate the restrictions 
listed directly after equation (1). Since the only equations which occur herein arise as 
a result of applying transformations of these types to equations of the form (29) to (31), 
these restrictions will always be met . 

It will be proved subsequently that a change of dependent variable of the form (26) 
applied to any equation in one of the canonical forms (29) to (31) transforms this equation 
into one which has the same canonical form (only the coefficients a, b, and c are 
changed). Since the order in which the transformations (26) and (27) are applied to a 
given equation is immaterial, the combined effect of applying a transformation of the 
type (26) and one of the type (27) to a given equation can be analyzed as follows. First, 
determine what conditions must be satisfied by the coefficients of an equation which is in 
one of the canonical forms (29) to (31) and which in addition can be transformed into a 
weakly separable or a separable equation by a change of variable of the type (27). Once 
this is done, it is only necessary to decide which of the differential equations having this 
canonical form can be transformed by a change of variable of the type (26) into an equa- 
tion whose coefficients satisfy these conditions in order to determine which of the equa- 
tions in this canonical form can be transformed by a combined change of variable into a 
weakly separable or a separable equation. 
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TRANSFORMATIONS OF THE DEPENDENT VARIABLES 


Apropos of these remarks, we now turn to a discussion of how equations of the types 
(29) to (31) transform under a change of variable of the type (26). This material is en- 
tirely equivalent to that given in reference 3. However, it is convenient to rederive some 
of the results in a form which is more suitable for our purposes. 

First, consider the hyperbolic equation (29). Substituting 

V(x , y) = A(x,y)U(x,y) 


into this equation yields 


xx 


- V + aV 

yy a 


+Ev y + cv = o 


(32) 


where 


a = a - 1 A (33) 

X x 

b = b + — A (34) 

X ^ 



This proves that, when a transformation of the type (26) is applied to an equation of the 
type (29), the form of the equation is unaltered. Evidently, 
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y 







r^> 





Hence, 

c - I (a +b ) - I (a 2 - b 2 ) = c - i (a +b ) - I (a 2 - b 2 ) (36) 

2 x Y 4 2 xy 4 

and 


Now define 





+ U - a + b 
x y x 



+ b 


x 


4 = c - 1 <a * + V ' j <a2 ' b2) 


(37) 


(38) 

(39) 


Equations (36) and (39) then show that, unlike the coefficients themselves, the quantities 
S H and / H are unaltered when a tranformation of the type (26) is applied to an equation 
which has the form (29). They will therefore be called the canonical invariants for an 
equation of the hyperbolic type. 

Now consider the parabolic equation (30) . Substituting 

V(x,y) = A(x,y)U(x,y) 


into this equation yields 


V xx + +SV y + = 0 


(40) 


where 


a = a 



(41) 
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(42) 


b = b 


c = c 



' b x x y 




(43) 


This proves that when a transformation of the type (26) is applied to an equation of the 
type (30) the form of the equation is unaltered. Evidently, 


Hence , 


t** 1 ~2 1 ~ 

c--a--a=c 
4 2 x 


l, 2 .I s .b»x 

4 . 9 x > y 


I (2c - la 2 - 

- a = 

foe - — a 2 - a \| 

b V 2 7 

y 

X 

b V 2 x / 


- a 


-ix 


Therefore, it follows from equation (42) that 



Equations (42) and (44) show that the quantities Jp and Jp are unaltered when a trans- 
formation of the type (26) is applied to an equation which is in the form (30). They shall 
be called the canonical invariants for the parabolic equations. 

Finally, consider the elliptic equation (31). Substituting 
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V(x,y) = X(x,y)U(x,y) 


into this equation yields 


xx 


+ V + aV 

yy x 


+ EV y + cV = 0 


(47) 


where 


a = a - 1 A (48) 

X x 

“b = b - — X (49) 

X y 

s— 

This proves again that, when a transformation of the type (26) is applied to an equation 
of the type (31), the form of the equation is unaltered. 

As before, it again follows easily from equations (47) to (50) that the quantities 

= a y " b x ( 51 ) 

/E = C -^ a x + V-7( a2+b2 ) (52) 

2 J 4 

are unaltered when a transformation of the type (26) is applied to an equation which has 
the form (31). They will be called the canonical invariants for the equations of the elliptic 
type. 

We have therefore shown that 

(C6a) An equation which is in one of the canonical forms (29) to (31) is transformed 
into an equation of the sa me form by any change of dependent variable of the type (26). 

(C6b) For each of the canonical forms (29) to (31), there are two quantities, called 
the canonical invariant s, which are unaltered when a transformation of the type (26) is 
applied to an equation which has that canonical form. 


X / X N 

-L lb- 2 ^ 
A V X, 


(50) 


TRANSFORMATIONS OF THE INDEPENDENT VARIABLES 

Having discussed the effects of transformations of the type (26) on the canonical 
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forms (29) to (31), we now turn to a discussion of the effects of transformations of the 
type (27) on these canonical forms. To this end notice that the three canonical forms (29) 
to (31) can be combined into the single equation 


£U +i ^U +a 3U +b 3U + cU = 0 

for 

i = -i, o, l 

(53) 

9x 2 9y 2 9x ^ 




It is convenient to define the linear operator for 

j = - 

1, 0, 1 by 


L<1)(U) . a2u + j s2u + a su * b 3U 

for 

j = -i, 0, 1 

(54) 


9x 2 9y 2 dx dy 

When the general transformation 


€ = <p(*, y) 

i = ^(x,y) 


H<pju) * o 
3(x, y) 


(55) 


of the type (27) is applied to equation (d 3), an equation of the form (1) with discriminant 
2 

D = /3 - ay is obtained. The coefficients of this equation depend only on the coefficients 
of equation (53) and the functions cp and \p. In order to determine this dependence, it 
is only necessary to apply the chain rule and then substitute the results into equation (53). 
Thus, 


U x = Vx + Vx 

Uy = V y + U^ y 

U xx - U |^x + 2U ^x + Vl + + 

u yy = u l^y + 2U + u w^y + u ? <p yy + u ^yy 

Substituting these results into equation (53) and collecting terms reveals that the coeffi- 
cients in equation (1) are determined by the following equations: 

a =<p 2 + ]cp 2 (56) 

0 = <?Y^ X + j^y^y ( 57 ) 
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y = + My 

(58) 

A = L (j) (<p) 

(59) 

B = 

(60) 

C = c 

(61) 


TRANSFORMATIONS WHICH LEAD TO SEPARABLE EQUATIONS 

Equation (1) is separable if and only if its coefficients satisfy conditions ( 8 ) to (13). 
Hence, it follows from equations (56) to (61) that equation (53) can be transformed into a 
separable equation by a change of variable of the type (55) if, and only if, 


fd 1= ^ +j? ,2 (62) 

0 = ?x*x + i<p y^y ( 63 ) 

fej = + j^y (64) 

fd 2 = L^(<?) (65) 

fe 2 = L ( 66 ) 

f(d 3 +e 3 ) = c (67) 


where the d k for k = 1, 2, 3 are functions of £ = <p(x,y) only and the e k for k = 1, 

2, 3 are functions of 77 = i^(x,y) only . It has already been pointed out that the sign of the 
discriminant of the equation of type (1) which results from the transformation (55) must 
be the same as the sign of the discriminant of equation (53). Hence, it follows from con- 
ditions (15) to (17) that 


dj > 0 and e^ < 0 

if equation (53) is hyperbolic (i.e. , if j = -1), 

dj + 0 and e^ = 0 


( 68 ) 


(69) 
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if equation (53) is parabolic (i.e. , if j = 0), and 


dj > 0 and e x > 0 (70) 

if equation (53) is elliptic (i.e. , if j = 1). These conditions are, in fact, direct conse- 
quences of equations (62) to (64). 

On the other hand conclusion (C2) shows that equation (1) can be weakly separable 
but not separable only if it is hyperbolic. Since the type of equation cannot be changed 
by a transformation of the form (27), we conclude that equation (53) can be transformed 
into a weakly separable but not separable equation by a transformation of the type (27) 
only if it is hyperbolic. Now equation (1) is weakly separable but not separable if, and 
only if, its coefficients satisfy condition (25) and either conditions (19a) to (24a) or con- 
ditions (19b) to (24b). In view of the symmetry of these conditions, however, we can as- 
sume without loss of generality that only conditions (19a) to (24a) are relevant. Hence, 
it follows from equations (56) to (61) that equation (53) c an be transformed into an equa - 
tion which is weakly separable but not separable b y a change of variable of the type (55) 
if, and only if, it is hyperbolic and 


0 = cp 2 - cp 2 
^y 

(71) 

fd l = Vx " ^y * 0 

(72) 

, ,2 ,2 
fe l =x f y x ~ ^y 

(73) 

o = l( - 1 V) 

(74) 

f(e 2 +d 2 ) = L (_1) (^) 

(75) 

fe 3 = c 

(76) 


At this point, it is convenient to consider the three equations (29) to (31) individually. 
We discuss the hyperbolic equation (29) first. 


EQUATIONS OF HYPERBOLIC TYPE (j = -1) - SEPARABLE CASE 
Functional Form of Invariants 

First, suppose that equation (29) can be transformed into a separable equation by a 
change of independent variable of the type (55). Then the functions cp and \p must sat- 
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isfy conditions (62) to (67) with j = -1 and dj > 0 and < 0. It is therefore permis- 
sible to introduce the (real) functions 



and 

1 

and to define a (real) function u of $ only and a (real) function v of 77 only by 



Thus, there are functions u and v of x and y such that 

u(x,y) = u[^(x,y)] =u(£) 
v(x,y) = v|V(x,y)] = v(?y) 

We shall suppose that these equations can always be solved for <p and i// (if necessary 
the domain of the differential equations can be divided into a series of subdomains in 
which these equations can be solved) to obtain 

£ = <p(x,y) = ?[u(x,y)] (79) 

and 

V = ^(x, y) = ^[v(x, y)] (80) 

Therefore, it follows from equation (77) that 
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In a similar way, we find 


Hence, 


Similarly, 


i£ = VdT — 
ay T 1 ay 

(82) 

dif/ _ r— dv 

V e l^ 

(83) 

_ ^/li - 9v 

(84) 



(85) 




+ — d* 






9y 2 9y 2 2 W 


( 86 ) 

(87) 

( 88 ) 


It now follows from definition (54) and equations (81) to (88) that 



(89) 


L (_ 1 ) (<p) =-d\ 
2 1 


S-'k*) = 


/au\ 2 _ /au\ 2 
\3x ) ysyy 


.(-1)/ 




dv ' 
i3xj 


w 


V 1 ®! 


L (_1) (v) 


(90) 


Substituting equations (81) to (84), (89), and (90) into conditions (62) to (67) with j -1 
yields 


, 2 2 
f = u - u 
X y 


0 = Vx - u y v y 


-f = v 2 - v 2 
X y 


fd 2 d i ("2 - 4 ) + L 


(- 1 ) 


(U) 


^ 2 = -iei(v 2 -v 2 ) + 


(- 1 ) 


(v) 


(91) 

(92) 

(93) 

(94) 

(95) 


f(dg + e 3 ) = c 


(96) 


Upon substituting equation (91) into equations (94) and (96) and equation (93) into equa 
tion (95) we find 


L ( - 1, («)-d 4 (^-uj) 

!,<-%) = -e 4 (v 2 -v 2 ) 

c “< d 3 +e 3>( u x- u y) 

where the function d 4 of £, only and the function e 4 of r) only are defined by 


(97) 

(98) 

(99) 


d 4 = 


d 2 4 d i 


V 5 ! 
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In view of equations (79) and (80) we see that there are functions p 1 and p 2 of u only 
and functions and q 2 of v only such that 

Pl(u) = d 3 (|) = d 3 [£(u)] 
p 2 ( u ) = d 4 [^(u)] 

qi( y ) = e 3p^ v )] 

q 2 ( y ) = e 4 p(v)] 

Substituting these into equations (97) to (99) yields 

L ( " 1) (u) = (u| - u^p 2 ( u ) 

L ( ' 1) (v) = -(v^ - v|q 2 (v) 

c = [pj(u) + q 1 (v)]^u| - Uyj 
Equations (91) and (93) show that 


( 100 ) 

( 101 ) 

( 102 ) 


or multiplying both sides by v 


2 2 2 2 

U - U = V - V 

X y y x 


' (w) 2 ” v y( v y ' v x) 


Eliminating v y u y between this equation and equation (92) yields 

2 2 2 2 2 / 2 2 \ 

Vx - Vx = v y\y " v x) 


(103) 


or 
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2 ( 2 2 \ 2/2 2 \ 
M v y " v xj = v y\y - v xJ 


and since v y - v| = f * 0, this implies 



Hence, 

u* = ±v y (104) 

Substituting this result into equations (92) and (103) shows that 

v x = ±u y (105) 

where the plus sign in equation (105) must be associated with the plus sign in equa- 
tion (104), and the minus sign in equation (105) must be associated with the minus sign in 
equation (104). 

Differentiating equation (104) with respect to x, differentiating equation (105) with 
respect to y, and subtracting the results yields 

u xx " u yy = 0 < 106 ) 

Similarly, 


v - v =0 
xx yy 


(107) 


The most general solutions of these two equations are 


u = 1 F(x + y) + - G(x - y) 
2 2 


(108) 


v = - F(x + y) +1 G(x - y) 


(109) 


where F and F are any functions of x +y, and G and G are any functions of x - y. 
However, upon substituting equations (108) and (109) into equations (104) and (105) we 
find that 
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F* + G’ = ±(F' - G’) 


and 


F* - G’ = ±(F’ + G') 


or 


F» = ±F’ 


and 


-G’ = ±G' 


Hence , 


v = ± 


— F(x + y) - — G(x - 
2 2 


y) 


( 110 ) 


Since the general form of conditions (100) to (102) remains unaltered if we replace 
v by -v, no generality will be lost if we assume that only the plus sign holds in equa- 
tion (110). Hence, we conclude that the functions u and v must have the following 
form: 


U = i F(a) +1 G(t) 

(HD 

2 2 

v = I F(ct) - A G(t) 

(112) 


2 2 


where 


a = x + y (113) 

r=x-y (114) 

and F and G are any non constant functions of their arguments. (The functions F and 
G are nonconstant because f = (u^ - = F’G' #0.) 

It follows from definition (54) and equations (100), (101), (103), (106), and (107) that 
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au x +bu y = P2 (u) ( u x" u y) 

av x+ bv y = q 2 (v)(uj-u2) 

Multiplying the first of these equations by v and the second by u and then sub- 

J J 

tracting the resulting expressions yields 

a (Vy - v x u y) = (p2 v y - Vy) ' u y) 

By eliminating v y and u y from this equation by using equations (104) and (105) we find that 

a (4 - u y) = ( p 2 u x - ^2 V x)( U x " u y) 


or since 


, 2 2 , n 
f = Ux - u y * 0 


a = P 2 (u)u x - q 2 (v)v x (115) 

In a similar way, we find that 

b = -[p 2 (u)u y - q 2 (v)v y ] (116) 

Equations (102) and (111) to (116), in which pj, p 2 , q x , q 2 , F, and G can be any func- 
tions of their arguments, now give the most general form that the coefficients of equa- 
tion (29) can have if this equation is to be transformable into a separable equation by a 
change of variable of the type (55). For the present purpose, however, it is more con- 
venient to work with the canonical invariants of equation (29) rather than with its coeffi- 
cients themselves. To this end we differentiate equations (115) and (116) with respect to 
x and y to obtain 


a x = P2 (u)u xx ' <*2< v ) v xx + P2 (u)u x - *2 (v)v x 
a y = P 2 < u Ky - *2< V >V + P 2 (u)u x u y ' ^Yy 
b x = -[p 2 (uK y - q 2 (v ) v xy +p 2 (u)u x u y “ ^ v K v y] 

b y = -[P 2 ( u ) u yy - q 2 ( V ) V yy + P 2 (u)U y " q 2 (v)V y] 
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These equations together with equations (103), (106), and (107) show that 


a y +b x = 0 (117) 

and 

a x+ b y = [ P 2 (U) + q 2 (v) ](4 ' u y) ( 118 ) 

Also equations (92), (103), (115), and (116) show that 

a 2 - b 2 , {[p 2 (u)] 2 - [q 2 (v)] 2 ](u| - u 2 ) (119) 

Substituting these results together with equation (102) into the definitions (38) and (39) 
of the canonical invariants of equation (29) shows that 

/ H = 0 (120) 


and 






Or upon defining the function $ of u and the function $ of v by 

$(u) = Pj(u) - I p£(u) - - [p 2 (u)J 2 
2 4 

*(v) = qj(v) - i q 2 (v) [q 2 (v)] 2 


( 121 ) 


equation (121) becomes 


/ H = [*(u) + *(v)](uj - nfj (122) 

Thus, equation (29) can be transformed (by changing its independent variables) into a 
separable equation only if there are functions $ and $ of u and v, respectively, such 
that its canonical invariants satisfy conditions (120) and (122). 

Now suppose that equation (29) can be transformed into a separable equation by 
changing both its dependent and independent variables. Since the order in which these 
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transformations are performed is immaterial, it must first be possible to transform this 
equation by a change of variable of the type (26) into an equation whose coefficients satis- 
fy the conditions (102), (115), and (116). But this implies that equation (29) can be trans- 
formed into an equation whose canonical invariants satisfy conditions (120) and (122). 
Since the canonical invariants of equation (26) are unaltered by a transformation of the 
type (26), we conclude that if equation (29) can be transformed into a separable equation 
by changing both its dependent and independent variables , then it is necessary that its 
canonical invariants satisfy conditions (120) and (122). 

In order to see that these conditions are also sufficient, suppose that there exist 
functions 4> and St such that conditions (120) and (122) hold with u and v given by 
equations (111) and (112) for some functions F and G. Then it follows from definition 
(38) that condition (120) implies that there exists a function <jo such that 

(123) 

and therefore definition (39) shows that 



It is now easy to see from equations (32) to (35) that in this case the change of variable 
(ref. 3) 



V = e“/ 2 U 

transforms equation (29) into the equation 

V xx- V yy +/ H V = 0 


(124) 


(125) 


or, substituting equation (122), 

V xx - V yy + [*( u ) + *( v )]( u x " U y) V = 0 (l 26 ) 

Upon introducing the new independent variables u and v defined by equations (111) to 
(114) we find that 


V - V = (V - V )F’G' 
v xx yy ' v uu V vv' 
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and 


u 2 - u 2 = F’G' 
x y 

Substituting these results into equation (126) shows (since 


V - V + 
uu vv 


[# (u) + ’F (v)] V 


F’ * 0 and G’ * 0) that 
0 


(127) 


(128) 


and this equation is certainly separable. 

If any functions cp and which satisfy equations (79) and (80) were used as the new 
independent variables in place of u and v, then instead of equation (128) we would have 
arrived at the equation 


1 _d_ 

u'(«) 3| 



V'(V) 


_ 3 _ 

d-q 


z — V 

\v'(n) v j 


j<J?[u(|)] +^[v(77)]j V = 0 


(129) 


where u and v are the solutions of equations (79) and (80), respectively, for u and v 
in terms of | and rj, respectively. It is easy to see that this equation is also separable. 
We have therefore established the following conclusions: 

(C7) The canonical hyperbolic differential equation (29) can be transformed into a 
separable equation by changing both the dependent and independent variables if, and only 
if, there exist functions <f> and ^ and nonconstant functions F and G such that the 
canonical invariants and satisfy the following conditions: 


where 


> H - 0 

/ H =[$(u) + *(v)](u 2 - u 2 ) 


(130) 


u =- F(cx) +1 G(t) 
2 2 





a = x + y and r = x - y 


(131) 
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(C 8 ) If the canonical invariants of equation (29) do satisfy conditions (130) and (131), 
then this equation can always be transformed into a separable equation by introducing 
both a new dependent variable defined by 


V = e" /2 U 

where w is determined to within an unimportant constant by 



and new independent variables g and 77 defined by 

£ = cp( u) 

7i = $(v) 

where <p and 4 / are any convenient nonconstant functions (i.e. , cp' +_ 0 and 0 ) and 
u and v are determined from condition (131). 


Direct Calcu lational Procedure for Testing 

In practice, it may not always be easy to tell simply by inspection whether the invar- 
iant /pj of a given equation can be put in the form (131). Also, there may be several 
ways in which a given function can be expressed in the form (131). Each of these 
ways will lead to a different "coordinate system" in which the equation is separable. 

For these reasons, it is useful to give an alternative form of condition (131) which sup- 
plies the means of testing the invariant by direct calculation and which, in addition, 
gives a procedure for calculating all of the functions u and v which determine the new 
independent variables. 

To this end suppose first that <^H satisfies condition (131). Then equation (127) 
implies that the first equation (131) can be put in the form 


F'G’ 


= $ (u) + (v) 


(132) 


There exist functions $ and ^ such that equation (132) holds if, and only if, 
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I 


/] 


H 


3u av VF’G'y 


= 0 


But it follows from the middle two equations (131) that 


(133) 


-L = I F'((r) (— +— ^ 
da 2 \3u dv/ 

and 

± = .1 G '(t)(-L-±) 

dr 2 \3u 0v / 

Hence , 

a _ i _a_ + _i_ _a_ 

3u F’(cr) da G’(t) dr 

_d_ l _a_ _ l _a_ 

av F’(ct) da G'( r ) 3 t 

Using these in equation (133) shows that 


l a 

~i a //h| 

1 

a 

i a /^h\ 

F’ da 

F’ da\ F’GV 

G' 

dT 

G’ dr If'gv 


Hence, 


a 

j_ _a_(^i\ 

__d_ 

J_ a 

da 

F* acr \ F* / 

dr 

G* dT\G'/ 


Upon introducing the function S of a only and the function T of r only defined by 


S(a) = 



T(t) s 



(134) 
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this equation becomes 


2S 




+ 3T 'W t +T "/h 


(135) 


This equation was first introduced by Darboux in a less general context (ref. 3). 

Thus, if Sr satisfies condition (131), then there exist strictly positive functions S 
of ct only and T of r only such that Sr satisfies equation (135). Conversely, if there 
are strictly positive functions S of a and T of t such that Sr satisfies equa- 
tion (135), then it is an easy matter to reverse the steps carried out above to show that 
equation (134) can be used to define nonconstant functions F and G and equation (133) 
can be used to introduce functions <£ and such that Sr satisfies condition (131). 
Hence, we conclude that 

(C9) There exist nonconstant functions F and G and functions 4> and such that 
Sr satisfies condition (131) if, and only if, there exists a strictly positive function S of 
a only and a strictly positive functibn T of t only such that Sr satisfies equa- 
tion (1357 ! 

(CIO) If strictly positive function S and T ca n be found such that equation (135) 
holds, then the functions F and G for wh ich condition (131) is satisfied can be calcu- 
lated from equation (134) . 

Thus, if the functions S and T are known, conclusions (C8) and (CIO) give a pro- 
cedure for calculating a change in the independent variables which will transform equa- 
tion (29) into a separable equation. 

Now equation (135) always possesses at least one solution. (Note that equation (135) 
is satisfied by taking S = T = 0. ) We shall subsequently develop a procedure which will 
yield expressions for all the solutions S and T to equation (135) provided that Sr is 
not one of two special types of functions. These expressions will involve, at most, two 
undetermined constants. To determine what restrictions, if any, must be placed on these 
constants in order that these expressions satisfy equation (135), they must be substi- 
tuted back into that equation. If after this is done the constants can still be adjusted so 
that S and T are positive functions, then we can conclude that equation (29) can be 
transformed into a separable equation, and we can use the expressions for S and T to 
calculate the new independent variables which will accomplish this. 

Since this procedure will not work when Sr is one of two special types of functions, 
before establishing this procedure we shall prove that, if / H is one of these types of 
functions, the condition (131) is always satisfied and that the functions F and G can 
easily be determined. 

First, suppose that Sr = 0. Then it is easy to see that condition (131) can always 
be satisfied by taking 
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$ =V = 0 


and that any nonconstant functions F and G can be used for determining u and v. 
Next assume that + 0. We may then define the function in terms of by 



Br) 


for * 0 


( 1 - 


Suppose there is a constant c Q such that 




= c. 


(137) 


If c Q = 0 then 





But this implies that has the form 

/ H = 

for sone nonzero functions y and X. 

It is clear from equation (127) that condition (131) will always be satisfied if we take 


<f> = 1 


* = 1 


and 


F’ = y 
G’ = X 

If c Q + 0, definition (136) and equation (137) show that 
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But this is Liouville's equation. It is shown in reference 3 that the most general solu- 
tion to this equation is 


/„ = (138 

c 0 [?(o| - X(t)] 2 

for all non constant functions y and X. Now in view of equation (127), it is clear (if we 

o 

take y = F, X = G, and SI'(v) = -l/2c Q v ) that equation (138) satisfies condition (131). 
There are many other ways of choosing the functions F and G such that equation (138) 
can be written in the general form (131) and each of these choices, of course, leads to a 
different change of variable (coordinate system) which will transform equation (29) into 
a separable equation. To show this, notice that if y and X are bounded below we can 
always choose these functions in such a way that they are both positive, for in this case 
there exists a finite number M such that 


M = 


gib- {y(a),X(T)} 
cr, reD 


where D is the domain of definition of the differential equation. If we put 


y - ^ = M +y 
Xj = M + X 

then 

/ _ 2 

v — 

C 0 [y 1 (cr) - X^t )] 2 

and 


y x — 0 
XjSO 


A similar argument holds if both y and X are bounded above. Suppose that y and X 
are either both bounded above or both bounded below. Then no generality will be lost if 
we assume that y and X are both positive. Hence, if we define the nonconstant func- 
tions F and G by 
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1 


F(a) 


G(r) 



~\ 


V 4t3 - *2* 


dt 


g 3 



dt 


V 4t3 - S2 4 • g 3 


(139) 


for any suitable positive constants gg and gg, then it is shown in reference 12 (p. 438) 
that 


y(o) 



X(t) = 



where jp is the Weienstrass ^-function (ref. 12). Substituting these equations into 
equation (138) shows that 





and the results of example 1 in reference 12 (p. 456) whos that this can be written as 



F'G* 

2c o 


Ki F i G )-Ki F -H 


It is now clear from equation (127) that satisfies condition (131) with $ = ^?/2c o and 
In addition, the functions F and G can be calculated from equation (139). 
There is usually no difficulty in telling simply by inspection if a given function is a solu- 
tion to Liouville's equation. Hence, we have shown that 


35 



(Cll) The condition (121) can always be satisfied if either the invariant or the 
function defined by equation (136) is equal to a constant. 

Now suppose that /h + constant, # constant and that the functions S of a only 
and T of r only are any two simultaneous solutions of equation (135). Then dividing 
equation (135) through by /h and applying the differential operator 


8 2 

da St 


- V’h^h 


to both sides of the resulting expression, we find (after some manipulation, which is car- 
ried out in appendix A) that S and T must satisfy the equation 


3k 1 3k 9 

2 — is+5k 1 S' = 2 — - T + 5k«T' (140) 

3a 1 3t * 


where we have put 


■ 4-^4 

k = / 4 9 / 

Notice that equation (140) can also be written more compactly as 



(141) 


(142) 


(If in any equation, functions which are negative are raised to fractional powers, 
it will always be possible to eliminate the fractional exponents by carrying out the 
indicated differentiations to obtain an equation in which all terms are real. There- 
fore, no difficulty will be encountered when this more compact notation is used. ) 

Equation (140) can be used to obtain an equation which involves T only and not S 
and an equation which involves S only and not T. Actually, this can be done in several 
ways. We shall discuss one of these ways in detail and also indicate briefly how one of 
the alternative methods can be carried out. 

If kg were zero, definition (141) would imply (since by hypotheses # 0) that 
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d A , 0 

3t 

or equivalently, that there exists a function y of a only (which by hypotheses is non- 
constant) such that 



Therefore, equation (141) would show that 

k 2 = ^H y ’ 


Substituting this expression together with the condition kg = 0 into equation (140), how- 
ever, yields (after some manipulation) 


4 9/ h V^hI 

2S — = 2S 

S H 3a *4 

* H 


2y M S + 5y^ 

y' 


If S * 0, then 


8 8 /h _ y* 'S + Sy'S 1 

/ H 3a y'S 


Now the right side is a function of a only; hence, this expression implies that 




= 0 


which is contrary to hypothesis. We therefore conclude that = 0 implies S = 0. A 
similar argument shows that k^ = 0 implies T = 0. 

We have therefore established the following conclusion: 

(C12) Whenever /h * constant and + constant, (/jj) = 0 implies that equa- 
tion (131) is satisfied only if S = 0, and = 0 implies that equation (131) is satisfied 

only if T = 0. 
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Since if S = 0 there is no need to calculate T and the reverse , we shall now make 
the additional assumption that kj # 0 and k 2 * 0. Since k 1 * 0, we can divide it into 
both sides of equation (140). Differentiating the result with respect to t yields 



= 0, equation (143) is essentially a second-order ordinary differen- 

T 

tial equation for T. (Notice, however, that the coefficients in this equation are, in gen- 
eral, functions of the two variables a and r and not just of r as is ordinarily the 

=* 0, we can divide it into both sides of equation (143) to obtain 

T 


case.) If l/kjfkj 




Equation (144) can now be substituted into equation (142) to obtain, in this case also, what 
is essentially a second-order linear ordinary differential equation for T. In either case 
then, it follows from equations (142) to (144) that if we define t^, t^, and t 2 ' by 
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then 


'T +t 


(^T* + t^T" = 0 


(148) 


In view of the symmetry of equation (140), it is easy to see by interchanging a and 
t that if we define s^, s^, and s^ by 
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for Kg * 0 


sd> = j 


5 k 3/5 _9_ 
2 2 3 t 


K 2 \ k l, 


2/5 


AAk 2 / 5 
lr 1 

V 2/ 


- k. 


-Jaj 


5 1 j9_ 

2 ,2/5 9a 
K 1 


lk 2 ' 


for Kg = 0 


8® = < 


5 .3/5 3 

— Ko 


0T 


k l / 5 / k i' 


K 2 \ k 2, 


for Kg ^ 0 


— for K 9 = 0 
k * 

k 2 


with 



(150) 


(151) 


then 


s^S +s^^S* + s^S” = 0 (152) 

There is another way in which equations of the form (148) and (152) can be derived. 
This procedure will lead to equations whose coefficients t^ and s^ (i = 0, 1, 2) are, 
in general, different from those obtained above. To this end, equation (140) is differen- 
tiated first with respect to a and then with respect to t to obtain the following two 
equations : 


and 



+ SkjS" 


= 2 



T + 5 



(153) 




= 2(k 2 ) TT T + 7(k 2 ) T 


T’ + 


5kgT M 


(154) 
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Since kj * 0, kg * 0, and i= 0, in none of the equations (135), (140), (153), and 
(154) can the coefficients of T and its derivatives and S and its derivatives all be equal 
to zero. Hence, S, S', and S" can be eliminated among these equations. This proce- 
dure will result in an equation for T which has the same form as equation (148). Sim- 
ilarly, T, T', and T" can be eliminated between equations (135), (140), (153), and (154) 
to obtain an equation for S which has the form (152). This procedure has the advantage 
over the preceding one that one less differentiation is required to calculate the coeffi- 
cients of the equations. However, the algebra involved is more tedious. 

We have now proven that if is not of the form covered by conclusions (Cll) and 
(C12), then all simultaneous solutions T and S of equation (135) must satisfy equa- 
tions (148) and (152), respectively. We shall show how these equations can be used to 
determine S and T provided that the coefficients of equation (148) and the coefficients 
of equation (152) do not all vanish. 

In appendix B, it is shown that the following four statements are equivalent: 

(1) All the coefficients of equation (148) vanish. 

(2) All the coefficients of equation (152) vanish. 

(3) There exist functions X 2 and Ag of r only, functions y 2 and y 3 of a only, 


and a constant 


17 1 such that 




A^(t)/ 3 (ct) 


7 2 ( ct )] 5 [ [t 3 ( ct ) " * 3 ( t )] ' 


(155) 


and 


k 2 _ _7r l 


X 2 ( t ) 5 f X’ 3 (T)y^(a) 
y 2 ( CT )] 3 \ [^3 (a) “ X 3< t) _ 


(156) 


(4) There exist a function X 2 °f T only and a function y 2 of a only such that 


— = 
k l 


K 5 - 1 y 2 V 2 
K i o k i 

1 X 2 - 


(157) 
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Combining equations (155) and (157) and putting 


K = K 1 = K 2 


shows that 


K = -5 


Xjj( T )y^(CT) 

r 3 (a) - X 3 (t) 2 


Thus, K is a solution of Liouville's equation 



(158) 


(159) 


(160) 


It is shown in appendix C that in this case equation (142) always possesses nonzero solu- 
tions S and T which are given in terms of the functions appearing in equations (155) 
and (156) by 


Syo 7 t 2 2 

— = ^ r 3 +7r 3 r 3 +7r 4 

y o « 


(161) 


"2 2 

-TAoX’o = X„ +7ToXo + 7 7 




3" 3 


(162) 


where 77 2 to 77^ are arbitrary constants. 

Suppose, therefore, that t^, t^, and t^ are not all zero and, hence, that s^), 
s^, and s< 2 > are not all zero. If the ratio of each pair of nonzero coefficients of equa- 
tion (152) is a function a only, then, after division by a suitable factor, equation (152) 
becomes just an ordinary differential equation (its coefficients are functions of a only) 
which can always, in principle, be integrated to obtain an expression for S as a function 
of cr only. If this is not the case, then equation (152) can be divided through by a non- 
zero coefficient to obtain an equation which has the same form as equation (152) but which 
has the properties that one of its coefficients is equal to unity and at least one of its 
remaining coefficients is not independent of r. This equation can then be differentiated 
with respect to t to obtain an equation of the form 
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m 4* n and m , n = 0 , 1 or 2 


(163) 


d m S 


+ s 9 ^§ = 0; 


da 


m 


do» 


and 


* 0 (164) 

If Sg/s j is a function of a only, then this equation can be solved to determine S as a 
function of a only. (If m = 0 and Sg = 0, equation (163) would then show that equation 
(135) possesses only the solution S' =0.) But if Sg/ij depends on t, then this equa- 
tion shows that 


^nin jIIq 

= = 0; m * n and m, n = 0, 1 or 2 (165) 

da m da 11 

If in this case the smaller of the two integers m and n is zero, we conclude that equa- 
tion (135) has only the solution S = 0. 

We have therefore shown that, whenever is not of the form covered by conclu- 
sions (Cll) and (C12) or by equations (155) and (156), all the solutions S to equation (135) 
can be determined to within at most two arbitrary constants by solving the appropriate 
one of the three ordinary differential equations (152), (163), or (165). 

Similar considerations, of course, apply to the function T and equation (148). Thus, 
whenever /h is not covered by conclusions (Cll) and (C12) or by equations (155) and 
(156) all the solutions to equation (135) can be obtained by solving ordinary differential 
equations which are, at most, of sec ond order . 

This completes our discussion of the conditions for transforming equation (29) into a 
separable equation. 

EQUATIONS OF HYPERBOLIC TYPE (j * -1) - WEAKLY SEPARABLE 

BUT NOT SEPARABLE 

Functional Form of Invariants 

Now suppose that equation (29) can be transformed into an equation which is weakly 
separable but not separable by a change of independent variable of the form (55). Then 
the functions <p and i// must satisfy conditions (71) to (73). Differentiating equation (71) 
with respect to x shows that 
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Ill II III I III 


^xx^x ” ^xy^y ® 


(166) 


and differentiating equation (71) with respect to y shows that 


^x^xy ~ Vy^yy 0 


(167) 


Multiplying the first of these equations by <p x , the second by </> y and then adding 
the results shows that 


Vxx 


CD CD =0 

Vyy 


(168) 


or, using equation (71), 



But equation (72) shows that 


(169) 


cp + 0 and 

X 


Vy *0 


(170) 


We therefore conclude that 

^xx ~ ‘Pyy = 0 ( 171 ) 

Definition (54) and equation (74) now show that 

a<p x +bcp y = 0 (172) 

Equation (72) shows that d ^ + 0. We can therefore define a nonconstant function u of | 
only by 



Thus , there exists a function u of x and y such that 

u(x,y) = u[<p(x,y)] = u(£) 


(173) 


(174) 
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Hence , 


u„ 



and 


u 


y 



Using these results in equations (71) and (72) shows that 



(175) 


and 


f = V// x - U y i// y * 0 (176) 

2 

Substituting equation (176) into equation (73) and then adding -ej/4 times equation (175) 
shows that 



or, upon collecting terms, 



(177) 


Differentiating this equation with respect to x gives 




flu 

2 x y 



= 0 


and differentiating with respect to y shows that 
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= 0 


Multiplying the first of these by [^ x - (e^/2)] and the second by [^ y - (e 1 u y /2)J and 
then adding the results shows that 



Upon using equation (177) this becomes 



e'-, 


~ „ (u xx u yy^ 0 ^x u x " ^y u y^ 


yy 


y y' 


+ 




( \b u 
vv x y 


“x^ = 0 


Hence, upon using equation (175) and adding and subtracting (e 1 /2)u x and (e 1 /2)u y we 
find 



Equation (177) now shows that 



'/ \ 2 

/ 0 x 2 ! 


1 

2 


+ (*W“ y )_ 

^xx - ^yy - e'lWVx - VV 


Since equation (72) would be violated if both 
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and 




e 1 

\b i u =0 

2 y 

we conclude that 

+**-*„- e i<Vx - Vy> = 0 <178) 

Definition (54) shows that, when equations (176) and (178) are substituted into equation 
(75), we obtain 

(d 2 + e 2 - e’ 1 )(i// x u x - i/yi y ) = a^ x +bi// y (179) 

Equation (172) and definition (174) show that 

au x + bu y = 0 (180) 

Multiplying equation (179) first by u and then by u shows, after using equation (180), 
that x y 

(d 2 + e 2 - e’ 1 )u x (^ x u x - ^ y u y ) = b(u x ^ y - i^u y ) (181) 

and 

(d 2 + e 2 - e* 1 )u y (^ x u x - ^ y u y ) = a(<// x u y - <^ y u x ) (182) 

It follows from equation (175) that 

- Vy } = U y (l// x u y - Vx } 

and 

u y(^x u x - V9 = Wy " Vx ) 
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and equation (176) implies that 


xly n - \L/ u ^ 0 
y 7 x 

Substituting these results into equation (181) and (182) shows that 

b = ” U y( e 2 " e l + d 2 ) (183) 

a = \ (e 2 " e l + d 2 } ( 184 ) 

Finally, substituting equation (176) into equation (76) shows that 

c = ' u yV e 3 (185) 

Equations (183) to (185) with u and determined from equations (175) and (177), 
respectively, now give the most general forms that the coefficients of equation (29) can 
have if this equation is to be transformable into an equation which is weakly separable 
but not separable by a change of variable of the type (55). As before, we shall use these 
relations to obtain expressions for the canonical invariants. To this end, we differenti- 
ate equations (183) and (184) with respect to x and y to obtain 

a x = »W d 2 +e 2 " e l> + u x*x< e 2 “ e i)' +u x d 2 d l 

a y = u xy< d 2 +e 2 ‘ e 'l) + W (e 2 ' ^ +u x u y d 2 d l 

b x = -u xy (d 2 + e 2 - e ’l) - u y^ e 2 ' e V' ’ u y u x d 2 d l 

b y = " u yy (d 2 +e 2 " e l ) " u y^ / y (e 2 " e l)’ " u J d 2 d l 

Hence , 

a y +b x = K*y ~ u y^x^ e 2 " e i)' 

and since the same derivation which was used to derive equation (171) from equation (71) 
suffices to show that equation (175) implies 
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we see from equation (175) that 


u 


xx 


u 


yy 


= 0 


a x +b y = Mx - u y V (e 2 “ e i ) ’ 


Finally, equations (175), (183), and (184) show that 


a 2 - b 2 = 0 


Using these results together with equation (185) in definitions (38) and (39) shows that 

y H = " u y^x )(e 2 " e l ) ’ (186) 


/h = (u x^x - u yV 


e 3 ' \ (e 2 - e i)’ 


(187) 


Upon multiplying equation (177) through by (u^ and using equation (175) we find 

2 / \ 2 


u %- : 7 u *) ' U ?fy-T'v) =0 


Hence, upon factoring and using equation (175) again 


0 = 


M*x- T u x +u ylVT U y 


“xl^x-y'Sc +u ylVT u yj 


K^X - U y^y) 


Therefore, equation (176) shows that 


M^x - ~ u x) +u y K 



(188) 


Let e 4 be any convenient nonzero function of r] and define the nonzero f miction e + of 
7] by 
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for # 0 


(189) 


e + = 


for e j = 0 


It is also convenient to define two new functions eg and eg of rj only by 

e 5 = e+ ^ e 2 “ e i )f 


” e+ If 3 " 2 (eg ~ ^ 


Then since e + # 0, equations (186) and (187) become 


J H = K*y ~ u y+x> — 


S H = - u yV 


And again since e + + 0, we can define a nonconstant function v of rj only by 


J e + (r 


v(v) = 2 / — — drj 

e + (rj) 


This shows that there are functions v and w of x and y such that 

v(x , y) = v[i^(x , y)] = v(r/) 
fv(x,y) - u(x,y) for e x * 0 

v(x,y) for ej = 0 


w(x,y) =■ 


Then 


V - v'l// = xp 
X Y x r x 


(190) 

(191) 

(192) 

(193) 
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And definitions (189) and (193) show that 


w„ 


7 w y 



= \L> i u 

2 y 


Substituting these results into equations (188), (190), and (191) shows, after using equa- 
tion (175) and the fact that e + * 0, that 

U x w x +Uy w y = ° (195) 

y K=\ (u x w y ' u y w x )e 5 < 196 ) 

=\ K w x ~ u y w y )e 6 (197) 

Equations (175) and (195) now show that either 

Ujj = u y and w x = -w y 


= -Uy and w x = w y 

Hence, if we put 

a = x + y'' 

► 

t = x - y> 


(198) 


then we conclude that there exists a nonconstant function F of cr only and a nonconstant 
function G of r only such that 
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u(x,y) = F(cr) and w(x,y) = G(r) 


or 


(199) 


u(x,y) = G(t) and w(x,y) = F(a) 

Upon substituting these results into equations (196) and (197) we find that 


✓h - -F’G'e^ 




V for u = F and w = G 


/ H = F'G’e, 


y H = F'G'e 5 


for u = G and w = F 


/ H = F ’ G ’ e 


6 


J 


Now equation (192) can be solved for r] as a function of v to obtain 

V = ?(v) 

and equation (174) can be solved for £ as a function of u to obtain 

£ = $>(u) 

We can therefore define functions <f> and ^ of v only by 


$(v) = 


- e 5 (??) = -e 5 | < P( v )J for u = F and w = G 




*(v) s e e (rj) = e 6 p(v)] 


for 

u - G 

and 

II 

£ 

r for 

u = F 

and 

O 

II 

£ 

I for 

u = G 

and 

O 

II 

£ 


( 200 ) 


( 201 ) 


( 202 ) 


Substituting these into equation (200) shows that 
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J H = $(v)F'G’ 
/ H = *(v)F'G’ 


(203) 


Equations (193) and (199) show that v must be given by one of the three following 
equations: 


and that 


when equation (204a) holds. 


when equation (204b) holds, and 


v = F(cr) + G(t) 

(204a) 

v = F(ct) 

(204b) 

v = G(t) 

(204c) 

= F or u = G 

(205a) 

u = G 

(205b) 

u = F 

(204c) 


when equation (204c) holds. The new independent variables are then given by equations 
(201) and (202). 

Thus, equation (29) can be transformed into an equation which is weakly separable 
but not separable (by changing the independent variables) only if there are functions 4? 
and ^ of v such that the canonical invariants satisfy equations (203). Now the same 
argument that was used in the separable case suffices to show that if equation (29) can be 
transformed into an equation which is weakly separable but not separable by changing 
both its dependent and independent variables then it is necessary that its canonical invar- 
iants satisfy conditions (203). 

In order to see that these conditions are also sufficient suppose that there exist func- 
tions 4> and such that conditions (203) hold with v given either by equation (204a) 
or (204b) . Then it is easy to verify by direct calculation that 
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- - I G'(r) f $ (v) dvl = — --+-G'(r) f 4>(v) dv 
*» - 1 J 9x L 2 4 •* 


0y L2 4 


(i) 


Hence, there exists a function oj v ’ such that 


^(i) _ a _ G t( r ) r $( v ) dv 

x 2 J 

= -- +- G'(t) f $(v) dv 

y 2 4 J 


for v = F + G and v = F 


(206) 


On the other hand, if conditions (203) hold with v given either by equation (204a) or 
(204c), then it can again be verified by direct calculation that 


_a_ 

3y 


— +-i F'(ct) J* $(v) dv =— + |f'(o) J' 4>(v) dv 


d_ 

9x 


2 4 


( 2 ) 


Hence, there exists a function co v ' such that 


=— +— F’(o') f $(v) dv 

x 2 4 ^ 

w (2) _ _b + _1 J" $( v ) dv 


> for v = F + G and v = G 


(207) 


It now follows from equations (32) to (35) that in this case the change in variable 


JX) 

V = e w U 


(208) 


transforms equation (29) into the equation 


\ G’[/ dv](V x + V y ) + (> H *1 y H )v = 0 


V - V +• 
xx yy 


when v = F + G and when v = F and that the change in variable 


w (2) 

V = e w U 


(209a) 


transforms equation (29) into the equation 
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(209b) 


V - V P- [/ * (V) dv](V x - V y ) + (/ H - | •'h) V ’ 0 
when v = F + G and when v = G. Hence, upon substituting equations (203) we obtain 

V xx " V yy + \ G ' [y dv ]( V x + V + F ’ G ' (* + | V = 0 (210a) 

for v = F + G and v = F, and 

V xx " V yy ' ^ F ’ [y *(▼) dy ](V x - V y ) + F’G ’ - i ^ V = 0 (210b) 

for v = F + G and v = G. Upon introducing the new independent variables v = F + G 
and u = G we find that 


V + V = 2F'V 
x y v 


V -V = 4F'G , V +4F'G'V 
xx yy uv vv 


and upon introducing the new independent variables v = F and u = G we find that 

V +V = 2F’V 
x y v 


V - V = 4F’G’V 
xx yy uv 


Substituting these results into equation (210a) shows that 


V uv +v vv + |[/* (v > dv ] V v + J 


't(v) + — 4>(v) V = 0 for v = F + G and u = G (211a) 
2 


uv 


+ i 


4 



$(v) dv 



+ 1 $(v)lv = 0 

2 


for v = F and u = G. 


(211b) 


Upon introducing the new independent variables v = F + G and u = F we find that 
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V - V = 2G T V 
x y v 


V -V = 4F'G'V +4F'G'V 
xx yy uv w 


and upon introducing the new independent variables v = G and u = F we find that 


V - V = 2G’V 
x y v 


V - V = 4F , G , V 
xx yy V uv 


Substituting these results into equation (210b) shows that 


V + V 
uv w 


-;[/ 


4>(v) dv 


V v + 1! 


ptf' (v) - - (v) 

2 


V = 0 for v = F + G and u = F (211c) 


and 


V. 


uv 


;[/ 


4>(v) dv 




*(v) -i$(v) 
2 


V = 0 for v = G and u = F 


(2 lid) 


It is clear that equations (211a) to (2 lid) are weakly separable but not separable. It is 
not hard to show that if, instead of taking u and v as the new independent variables, 
we had chosen any nonconstant functions of u and v as the new independent variables, 
we would have again transformed equation (29) into an equation which is weakly separ- 
able but not separable. 

We have therefore established the following conclusions: 

(C13) The canonical hyperbolic differential equation (29) can be transformed into an 
equation which is weakly separable but not separable by changing both the dependent and 
independent variables if and only if there exist functions $ and ^ and nonconstant func- 
tions F and G such that the canonical invariants and satisfy the following 
conditions : 
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>jj = <&(v)F , (o)G 1 ( T ) 

«/ H =^(v)F , (o')G , (t) 


where 


v = F(a) + G(t) 


or 

V (212) 

v = F(ct) 
or 

v = G(r) 

and 

a = x+y, T=x-y > 

(C14) If the canonical invariants of equation (29) do satisfy condition (212), then this 
equation can always be transformed into an equation whic h is weakly separable but not 
separable by introducing both the dependent variable V defined either by 


V = 



(213a) 


when v + F + G or when v = F or by 


V = 



(213b) 


when v = F+ G or v+G, where and co^) a re determined to within an unimpor- 

tant constant by equations (206) and (207), respectiv ely, and also the new independent 
variables | and rj defined by 


I = <p(u) 
V = $(v) 
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where cp and \p are a ny convenient nonconstant functions, v is given in equation (212^, 
u = F when the transforma tion (213a) is used and u = G when the transformation (213b) 
is used. 


Direct Calculational Procedure for Testing Invariants 

We shall now give an alternate form of condition (212) which can be used to test the 
invariants and by direct calculation. 

To this end notice that and satisfy condition (212) if, and only if, there 
exists a function F of a only and a function G of t only such that 



It is now an easy matter to establish the following conclusions: 

(Cl 5) There exist non constant functions F and G and functions $ and ^ such 
that and satisfy condition (212) if, and only if, there exists a nonzero function 
T of r only and/or a nonzero function S of a only suc h tha t either 


and 


9 (S/ H ) =±(T>„) 


da 


3t 


H' 






(214a) 


or 
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— (S> H ) = 0 

da n 


(S/ H ) = 0 


(214b) 


(T>„) = 0 


(214c) 


t <Vh> - 0 


(C16) If nonzero functions S and T can be found such that equation (214a) is satis- 
fied, then condition (212) holds with v = F + G, and the functions F and G are given by 


F' = — 
S 


(215a) 


G* = — 


(215b) 


If a nonzero function S can be found such that equation (214b) is satisfied, then condition 
(212) holds with v = G, and the function F is given by equation (215a). If a nonzero func- 
tion T can be found such that equation (214c) is satisfied, then condition (212) holds with 
v = F, and the function G is gi ven by equation (215b). 

It is clear that, if 4 0 and 4 0, then equation (214b) holds if, and only if, 


Ma ( /h )c 


= Function of a only 
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Similar remarks apply to equation (214c). Hence, there is no difficulty in determin- 
ing whether conditions (214b) and (214c) are satisfied. Determining whether and 
satisfy condition (214a) is slightly more difficult. If were zero, then it would 
be possible to transform equation (29) into a separable equation whenever it is possible 
to transform it into a weakly separable equation (compare condition (212) with conditions 
(130) and (131) and eq. (127)). Since this case has already been discussed, we shall 

were zero, then there would 


exclude it here and suppose that ± 0. 


be a function y of a only and a function X of r only such that 


S H = y(tf)Mr) 


This would imply that either 


T 


n l 


X 


or 


S 


*2 

y 


where irj and ^ are constants and that satisfies all three of the conditions (214a) 
to (214c). In this case then, it is only necessary to determine whether Sh satisfies any 
one of the following conditions in order to establish that condition (212) holds: 



Hence, suppose now that equation (214a) has a nonzero solution and that * 0 and 
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obtain 


We can therefore divide the first equation (214a) through by and 


S 



+ S* 



+ T' 


Upon applying the operator 



to both sides of this equation we obtain 

s H = t ('h) t 


(216) 


where 



d_ 

3 T 



Now by hypotheses ¥= 0. If there were a nonzero constant c q such that 


H 


= c. 


(217) 


then would satisfy Liouville’s equation and would therefore have to have the form of 
equation (138). Now it is clear (if we take y = F, X = -G, and <3>(v) = l/2c Q v^) that equa- 
tion (138) satisfies condition (212). Hence, assume that * constant. Thus, 

and do not both vanish. If one of them vanished, say (^jj) for definiteness, then 

equation (216) would show that T = 0; that is, that the first equation (214a) does not have 
a nonzero solution, which is contrary to hypothesis. We therefore conclude that 
(*h) * 0 an< f (*h) * 0* Equation (216) now shows (upon differentiation) that 
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11 

T 


(218) 


and 


S 


Nr 3 

(H 

Wa 5T 

P 

1 

Nc 3 

Nr 

(< h ) t 3° . 

hJ 


-n< 2 > 


= n^ 1 ) 


(219) 


These equations show that if equation (214a) has a nonzero solution, then must 

be a function of a only and must be a function of r only. Conversely, if these 

conditions are satisfied, equations (218) and (219) can be integrated to obtain 


S = 77^ 


da 


( 220 ) 


and 


f dr 

T = 7T 2 e (221) 

where 7 ^ and are constants, and these results can be substituted back into the two 
equations (214a) in order to determine whether they possess a nonzero solution. 

This completes the discussion of the hyperbolic equation (29). We shall now con- 
sider the parabolic equation (30). 


EQUATIONS OF THE PARABOLIC TYPE (j = 0) 

Functional Form of Invariants 

First, suppose that equation (30) can be transformed into a separable equation by a 
change of variable of the form (55). Then, the functions cp and if/ must satisfy conditions 
(62) to (67) with j = 0, dj + 0, and e^ = 0. In addition, we have already indicated that 
we shall require that b + 0 in order to avoid trivialities. Equation (64) now becomes 

if/ x = 0 (222) 
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(Incidently, this shows that equation (63) is automatically satisfied.) Equation (62) shows 
that 


fdj = <P*±0 

It follows from definition (54) and the fact that dj * 0 that substituting the two preceding 
equations into equation (66) yields 


b^y 


2 2 


<P 


x 


(223) 


Since by hypotheses 


d(<?>^) + o 
3(x,y) 

it follows from equation (222) that 

i/'y * 0 (224) 

Hence, equation (223) shows that 

e 2 * 0 (225) 

Equation (222) shows that i// is a function of y only and since 


0X 

we can conclude that e 2 is a function of 
zero function v of y only such that 


~z T x 

y only. It now follows that there exists a non- 


v 


e 2 




Hence, equation (223) becomes 
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(226) 


b 

v 


— <Px* 0 
d l X 


It can again be concluded from definition (54) and the fact that 
ting equation (223) into equations (65) and (67) yields 


dj * 0 that substitu- 


and 


^xx 


^2 2 

+ acp x+ bcp =—<p 

df 


c = (d„ + e«>) — cp 


2 

x 


(227) 


(228) 


Upon substituting equation (226) into equation (228) we obtain 

c_ d 2 +e 3 
b v 


(229) 


Since equation (226) shows that cp x * 0, we can divide equation (227) through by cp x to 
obtain 


a 



Equation (226) can be differentiated to obtain 



Upon eliminating ( P XX /<P X 


between these two equations we obtain 




(230) 


(231) 


Equations (226), (229), and (231), in which dj, d 2 , d 3 , e 3 , and v can be any function of 
their arguments, now give the most general form that the coefficients of equation (30) can 
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have if this equation is to be transformable into a separable equation by a change of vari- 
able of the type (55). However, it is again more suitable to use these expressions to 
obtain conditions for the canonical invariants. Since equation (45) shows that /p = b, 
we see from equation (226) that this is already partially accomplished. An expression 
for /p can be obtained by substituting equations (226), (229), and (231) into definition 
(46). After considerable algebraic manipulation, which is carried out in appendix D, we 
obtain 



(232) 


where the function R of cp is defined in appendix D. Definition (45) and equation (226) 
show that 


S p = b = JL <p\ * 0 
d. 


(233) 


This shows then that equation (30) can be transformed (by changing the independent vari- 
ables) into a separable equation only if there is a function cp oi x and y, functions dj 
and R of cp only, and a function v of y only such that the canonical invariants satisfy 
conditions (232) and (233). The argument used for the hyperbolic case now suffices to 
show that if equation (30) can be transformed into a separable equation by changing both 
its dependent and independent variables then it is necessary that its canonical invariants 
satisfy conditions (232) and (233). 

In order to see that these conditions are also sufficient suppose that there exist func- 
tions cp, R, dp and v such that equations (232) and (233) are satisfied, and define the 
function 8 by 


8 = J a dx 


(234) 


Then, 


a 


ae 

0X 


and definition (46) shows that 
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Thus, if we put 




a 2 - 






then definition (46) can be written as 



an 

dx 


Using this result in equation (232) gives 





r— — . 

a 

S2 - - R(ip) +- /— 1\ 

a 

h (—) 

dx 

v 2 Vv 

3 y 

_ \vj 


This shows that there exists a function co such that 


co 


X 



u v = Si - 1 R(cp) 
y v 2 \tp x ) 


On substituting equation (236) into equation (237) we obtain 

c 2 , 

S2 - w +— =- R{<p) 
y 2b v 

The results of appendix E show that equation (233) implies 


( 235 ) 


(236) 


(237) 


( 238 ) 
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_ i \ 


1 

2 b 

X 

i 


=1 b -Ibll 

2 y 2 V 

Substituting equation (236) into this expression shows that 


w =Ib -IbH 

“ 2 b x 2 y 2 V 


(239) 


It is shown in appendix F that equations (234), (235), (238), and (239) taken together 
with equations (40) to (43) imply that the change of variable 


v = e (l/ 2 )[u + (l/ 2 )ln|b| +®]u 
transforms equation (30) into the equation 


(240) 


+ - — )v + bV + -L 

R (a>) +- v' 

2 b J x y 2v 

2 


xx y 

Upon introducing the new independent variables 

£ = <p(x,y) 

V = y 

we obtain 


V = 0 


(241) 


+ 


( 1 b x^ 

^xx + h( Py - Vx K + 2 


V> +bV +JL 
S 71 2 v 


R(S) +-v'fa) 

2 


V = 0 


Hence, upon using equation (230) (which follows from differentiating eq. (233)) and equa- 
tion (236) we find 


9 dt 9 v 

<PZVtt +— cpZVt + bV + — 
x « 2 d, x ? 77 2 v 


RU) +- Vfa) 
2 


V= 0 


Finally, upon using equation (226) we obtain 
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V = 0 


+-d* 1 (|)V^ +v(r 1 )V ri +- 


R(0 + - v'fo) 
2 


and this equation is certainly separable . 

It is easy to see that equation (241) could also have been transformed into a separable 
equation if r/ was taken to be any nonconstant function of y. 

We have therefore established the following conclusions: 

(C17) The canonical parabolic differential equation (30) can be transformed into a 
separable equation by changing both the dependent and independent variables if and only 
if, there exist a nonconstant function cp, a nonzero function d^ of cp only, a function 
R of q> only, and a nonzero fun ction v or y only such that the canonical invariants 
J-p and Jp satisfy the following conditions: 



(242) 


(C18) If the canonical invaria nts of equation (30) satisfy condition ( 242), then this 
equation can always be transformed into a separable equation by introducing both a new 
dependent variable V defined by 


(1/2) [to +(l/2)ln|b| +0] 

V = e U 

where w is determined to within an unimportant constant by equations (236) and (237) 
and 6 is defined by equation (234), and new independent variables £ and ij defined by 

i = <p(x, y) 
n = ?(y) 

where cp is determined from condition (242) and $ is any nonconstant function of y 
only. 
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Direct Calculational Procedure for Testing / p 

We shall now give an alternate form of condition (242) which can be used to test the 
canonical invariants and To this end notice that (since <p x * 0) the second 

equation (242) can be written in the form 



Hence , upon differentiating by parts we can conclude that there exists a function R of 
cp only such that the second equation (242) holds if and only if, 



On the other hand, it is shown in appendix E that there exists a function d^ of cp only 
such that the first equation (242) holds if and only if, 


Upon putting 



1 

2 





in equations (243) and (244) we obtain 


(244) 
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Ml 


_a_ 

9y 


ii, T \ - v (; T )J = i [ T/p (l T \ - (; T ) y T ’ 


(245) 


and 


T +- JE T =- — (- 

x 2 b 2 b \v. 


(246) 


Since given T/v (recall that v * 0) the first-order linear differential equation 


<P x -^ y 


(247) 


always has a nonconstant solution (i. e. , a solution cp such that (p x and are not 
both zero) and since equation (247) shows that cp x = 0 and implies that <p y = 0, we can 
conclude that equation (247) always has a solution <p such that <p x ± 0, This remark 
taken together with the preceding results is sufficient to show that there exist functions 
ip , dp R, and v * 0 such that condition (242) is satisfied if, and only if, there exists a 
function T and a nonzero function v of y only such that equations (245) and (246) are 
satisfied. 

Multiplying equation (246) by (bT/v) shows that 


Hence, 



T (}L T 2 \ -(-t\ = (*) 

v \2v / x \v / y 2 \v/y 



(248) 


Using this result in equation (245) shows that 


d_ 

3y 





2 

and, by using equation (248) again to eliminate (bT /v) , we obtain upon differentiating 
by parts 
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( V H - K) y - ( T /p) x 


(249) 


The steps of this argument can easily be reversed to show that equation (249) taken 
together with equation (246) is equivalent to equation (245) taken together with equation 
(246). We have therefore established the following conclusions: 

(C 19) There exists a function cp and a nonzero function d^ of <p only, a function 
R of <p only, and a nonzero function v of y only such that Jp and Jp satisfy con- 
ditions (2 42) if, and only if, the reexist a function T and a no nzero function v of y 
only such that 

(yp) y • ( bT y) y ■ ( T *) x (249> 

1 b Y 1 b v 1 

T + ± _ T =- v-i - - v» (250) 

x 2 b 2 b 2 

(C20) If a function T and a nonze r o fun ction v of y only can be found such that 
equations (249) and (250) hold, then a function cp w hich s atisfies condition (242) can be 
found by solving the first-order linear partial diffe rential equation (247). Thus, cp is 
any function such that cp = Constant is an integral of the ordinary differential equation 

dx _ _ T 
dy v 

Notice that equations (249) and (250) are linear and homogeneous in T and v. 

Hence, they always possess the trivial solution T = v = 0. However, conclusion (C19) 
requires that the function v of y only be nonzero. We shall now develop a procedure 
which will yield expressions for all the solutions T and v to equations (249) and (250). 
As in the hyperbolic case, these expressions will involve undetermined constants. When- 
ever these expressions are obtained by this procedure, it is necessary to substitute them 
back into equations (249) and (250) to determine what restrictions, if any, must be placed 
on the undetermined constants in order that these equations be satisfied. If after this is 
done the constants can still be adjusted so that v # 0, then we can conclude that equa- 
tion (30) can be transformed into a separable equation, and we can use the expressions 
for T and v to calculate the new coordinates from equation (251). 

Now suppose that the function T and the function v of y only are any two simul- 
taneous solutions to equations (249) and (250). Then it is shown in appendix G that when 
the operator 
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d_ 

3x 


1 

2 



is applied to equation (250) the following equation is obtained: 

2v , / p+ v( /p ) y+ ^l = (^ p ) x T (251) 


where 



b / b \ ^ 

-_yy + Li ] (252) 

2b \2b/ 


Now if (jpj = 0> equation (251) is a third-order ordinary differential equation for 

v since this condition also implies that all the coefficients in this equation are functions 
of y only. On the other hand, if [jpj * 0, it can be divided into both sides of equa- 
tion (251) to obtain 


v 






(253) 


Equation (253) can be substituted into equation (250) to obtain, in this case also, what is 
essentially a third-order ordinary differential equation for v. (Notice, however, that 
the coefficients in this equation will be, in general, functions of the two variables x and 
y and not just of y.) In either case then, it follows from equations (250), (251), and 
(253) that if we define Y^, Y^\ and Y^ by 



(254) 
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i 


then 



r(0). 


v + Y^v' + Y' 


■(3)t 


= 0 


( 255 ) 


(256) 


(257) 


We have now proved that every function v of y only which satisfies equations (249) and 
(250) for any function T must satisfy equation (257). The coefficients of this equation 
cannot all be zero. For if (/p) = 0, this conclusion follows directly from equation (256); 

and, if {j-pj * 0, it is easy to see from equations (255) and (256) that Y'*^ = 0 implies 

that Y'*' = 5/2. Equation (257) can therefore be used to determine v in exactly the 
same way that equation (148) was used to determine T in the hyperbolic case. 

Thus, if the ratio of each pair of nonzero coefficients of equation (257) is a function 
of y only, then after division by a suitable factor equation (152) becomes just an ordin- 
ary differential equation (with coefficients now functions of y only) which can always be 
integrated to obtain an expression for v as a function of y only. On the other hand, if 
this is not the case, then equation (257) can be divided through by a nonzero coefficient to 
obtain an equation which has the same form as equation (257) but which has the properties 
that one of its coefficients is equal to unity and at least one of its remaining coefficients 
is not independent of x. This equation can then be differentiated with respect to x to 
obtain an equation of the form 


Yj^-^+Yg — = 0 m*n and m, n = 0, 1 or 3 (258) 

dy m dy 11 


and 
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Y 1 * 0 


If Yg/Yj is a function of y only, then this equation can be solved to determine v as a 
function y only. But if Y 2 /Y 1 depends on x, then equation (258) shows that 

H m .n 

= 0 m*n and m, n = 0, 1 or 3 (259) 

dy m dy n 

We have therefore shown that all the functions v of y which satisfy equations (249) 
and (250) can be determined to within at most three arbitrary constants by solving the 
appropriate one of the three ordinary differential equations (257) to (259). 

If (Jpj ~ equation (257) is just an ordinary differential equation (its coefficients 

are functions of y only). And, therefore, there are three linearly independent functions 
v of y only which satisfy this equation, and any linear combination of these solutions is 
also a solution. Now it is shown in appendix H that, for each of these infinitely many 
nonzero functions v of y only which satisfy equation (257), there exists a function T 
which contains two arbitrary constants, such that v and T satisfy equations (249) and 
(250). Thus, if (jpj = 0? there is a six-parameter family of solutions to equations 

(249) and (250). The analysis of appendix H shows that, in general, it may be necessary 
to use different expressions for T in different parts of the domain of equations (249) 
and (250). Since the analysis of appendix H is constructive, it can be used to determine 
T once the solutions v of equation (257) are found. As has already been shown, once 
T and v are determined the six-parameter family of independent variables which trans- 
form equation (30) into a separable can be found. 

If (/'pj * then once v is found by the procedure just described the function T 

can be determined from equation (253). It clear from the way that equation (257) was 
derived that, if v satisfies this equation, then T and v will automatically satisfy equa- 
tion (250). However, it is necessary to substitute these expressions for T and v back 
into equation (249) to determine what restrictions must be placed on the arbitrary con- 
stants which they contain in order that both equations (249) and (250) are simultaneously 
satisfied. 

Since a parabolic equation cannot be transformed into an equation which is weakly 
separable but not separable, this compleses the discussion of the parabolic equation (30). 
We shall now consider the elliptic equation (31). 
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EQUATIONS OF THE ELLIPTIC TYPE (j = 1) 
Functional Form of Invariants 


First, suppose that equation (31) can be transformed into a separable equation by a 
change of variable of the form (55). Then, the functions cp and ij/ must satisfy condi- 
tions (62) to (67) with j = 1, dj > 0, and ej > 0. It is therefore permissible to intro- 
duce the functions 


V 5 ! 


and 


1 

and to define a function u of £, only and a function v of rj only by 



(260) 


(261) 


Thus, there are functions u + and v of x and y such that 

u + (x,y) = u[<p(x,y)] =u(£) 
v(x,y) = v[i//(x,y)] = vfa) 

We shall suppose that these equations can always be solved for (p and ij/ to obtain 

I = <p(x, y) = ^[u + (x,y)] (262) 

and 

V = J//(x,y) = $ [v(x , y)] (2 63) 
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w~ 


Therefore 


Hence, 


It now 


, it follows from equations (260) and (261) that 


i£= V57— 

3x ’ 1 9x 

(264) 

i£= VdT— 

9y 3y 

(265) 

9x 1 dx 

(266) 

*t = */eT — 
3y ¥ 1 3y 

(267) 


= + i d , 

9x2 'ax 2 2 



£*.= WdTi^+Id'. 

9y2 9y 2 2 \9y 



9x 2 9x 2 2 W 


-2 , .2 . / . \2 

l±= jr*JL + ie\ l*L\ 

a 2 1 a 2 2 1 \ 3y / 

3y dy \ ■ y/ 


follows from definition (51) and equations (264) to (271) that 


L (1) (<?) =-dt 
2 1 




+ l(%+> 


L (1) (^) = — e i 

2 1 



+ V®I l( 1 V) 


(268) 

(269) 

(270) 

(271) 

(272) 

(273) 


76 


I 


Upon substituting equations (267) to (269), (272), and (273) into conditions (62) to (67) 
with j = 1 we obtain 


'-&)♦«) 

(274) 

0 = u x v x + U y V y 

(275) 


(276) 

[(<N(<f]**Ii < V> 

(277) 

fe 2 = | e l [W 2 + ( v y) 2 ] + V®1 l(1)(v) 

(278) 

f(dg + eg) = c 

(279) 

After substituting equation (274) into equations (277) and (279) and equation (276) into 

equation (278), we obtain 

L U) <»b = d 4 [«? + «? 

(280) 

L (1) (v) = e 4 [ ( v x) 2 + (v y ) 2 

(281) 

r 2 . . 2i 

c = ( d 3+ e 3 )[( u x ) +(u y y 

(282) 


where the function d 4 of | only and the function e 4 of r] only are defined by 
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Equations (274) and (276) show that 


«) + ( u y) * ( v x) 2 + ( v y) 2 
/ \2 

or multiplying both sides by (v y J 

Eliminating v^u^ between this equation and equation (275) yields 

(4) 2 [( v y) 2 + W 2 ] = ( v y) 2 [W 2 + ( v y) 

and, since { v x f‘ + (v y ) 2 = f * 0, this shows that 


(283) 


(0 - w ! 


Hence, 


^ = ±v . 


(284) 


Substituting this into equations (275) and (283) shows that 


u 


+ 

y 


= 


(285) 


where the minus sign in equation (285) must be associated with the plus sign in equa- 
tion (284), and the plus sign in equation (285) must be associated with the minus sign 
in equation (284). Thus, it is convenient to introduce a new function u of x and y by 


r \i + if the plus sign holds in equation (284) and 
the minus sign holds in equation (285) 

-u + if the minus sign holds in equation (284) and 
^ the plus sign holds in equation (285) 


( 286 ) 


Equations (284) and (285) now become 
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Equations (284) and (285) now become 


u = v 
x y 


u = - v 
y x 


(287) 


Equations (262) and (263) and definition (286) now show that we can define functions 
and P 2 of u only and functions q 4 and q 2 of v only such that 


Pj(u) = 


p 2 (u) s 


dg[<p(u)] if u = u + 
dg[^(-u)] if u = -u + 
'd 4 [^(u)] if u = u + 
d 4 [^(-u)] if u = -u + 
q 4 (v) = e 3 [^(v)] 
q 2 (v) = e 4 p(v)] 


Substituting these definitions together with definition (286) into equations (280) to (282) 
gives 


,^(u) = P 2 (u)(u 2 + U 2 ) 

(288) 

■ (1) (v) = q 2 (v)(v^ + v 2 ) 

(289) 

: [p 4 (u) +q 4 (v)](u^ +U y) 

(290) 


Now equations (287) are just the Cauchy-Riemann equations. Thus, u and v are 
conjugate harmonic functions. Therefore, 


V 2 U = V 2 v = 0 


(291) 


79 



and there exists an analytic function W of the complex variable 


z = x + ly 


( 292 ) 


such that 


W(z) - u(x,y) + iv(x, y) 


(293) 


In particular, we can conclude from this that 


dW 


dz 


= u 2 + u 2 = v 2 + v 2 * 0 
x y x y 


(294) 


It follows from definition (54) and equations (288) to (290), (291), and (294) that 


au x + bu y = p 2^ 


dW 


dz 


(295) 


av x + bv = q 2 (v) 


dW 


dz 


(296) 


c = [p 2 (u) + qj(v)] 


dW 


dz 


(297) 


It follows from the definition of the derivative of an analytic function that equations (295) 
and (296) can be written as the following single complex equation: 


(a + ib) — = (p 2 + iq 2 ) 
dz 


dW 


dz 


(298) 


Now equation (294) shows that dW/dz * 0. Therefore, it can be divided into both sides 
of equation (298) to obtain 


a + ib = (p 2 + iq 2 ) 


Hence, upon taking the real and imaginary parts we find 
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r 


= p 2 ( u ) u x + q 2 (v)v x 

(299) 

11 

CO 

'p' 

+ 

CO 

(300) 


upon recalling that the derivative of an analytic function is independent of direction. 

Equations (292), (293), (297), (299), and (300), in which pj, p 2 , q^, and q 2 can be 
any functions of their arguments and W can be any nonconstant analytic function of the 
complex variable z, now give the most general forms that the coefficients of equation (31) 
can have if this equation is to be transformable into a separable equation by a change of 
variable of the type (55). However, it is again more useful to use these expressions to 
obtain conditions for the canonical invariants. Upon differentiating equations (229) and 
(300) with respect to x and y we find, after using equations (291) and (294) to simplify 


a 


y 


- b„ 


= 0 


a x + b y = [P2< u ) + <l2 (v) 


dW 


dz 


Equations (287), (294), (299), and (300) also show that 


a 2 +b 2 =(p 2 + q 2 ) 


dW 

dz 


Substituting these results together with equations (297) into the definitions (51) and 
(52) of the canonical invariants of equation (31) shows that 


> E - 0 


and 




Or upon defining the function 4? 


of u and the function $ of v by 


(301) 


(302) 
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equation (302) becomes 


$(u) = Pj(u) - i p£(u) - I [p 2 (u)] 2 
^(v) = q x (v) - 1 qJ>(v) - i |q 2 (v)] 2 


/ E = [*(u) + *(v)] 


dW 

dz 


2 


(303) 


Thus, equation (31) can be transformed by a change in the independent variables into 
a separable equation only if there is an analytic function W of the complex variable z 
and functions $ and $ of u and v, respectively, such that its canonical invariants 
satisfy conditions (301) and (303). The argument used in the hyperbolic case now suffices 
to show that, if equation (31) can be transformed into a separable equation by changing 
both its dependent and independent variables, then it is necessary that its canonical in- 
variants satisfy conditions (301) and (303). 

In order to see that these conditions are also sufficient, suppose that there exists 
an analytic function W of the complex variable z and functions <f> and ^ such that 
equations (301) and (303) hold. Then it follows from definition (51) that condition (301) 
implies that there exists a function u> such that 


a = 



b = w y J 


(304) 


It is easy to see from equations (47) to (50) and equation (52) that the change of variable 
(refs. 3 and 4) 


V = e w / 2 U (305) 

transforms equation (31) into the equation 

V 2 V + / E V = 0 (306) 


where 


is the laplacian 
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Or substituting for in equation (303) this becomes 


V 2 V + [<f>(u) + 'lf(v)] 


dW 

dz 


2 

V = 0 


(307) 


But it is shown in reference 1 that, upon introducing the new independent variables u 
and v defined by equation (293), this equation transforms into the separable equation 

v uu + V vv + [* + * ( y )] V = 0 (307) 

It is also shown in reference 1 that, if any functions (p and i// which satisfy equations 
(79) and (80) were used as the new independent variables in place of u and v, then equa- 
tion (307) would still be transformed into a separable equation. 

We have therefore established the following conclusions: 

(C21) The cano nical elliptic differential equation (31) can be transformed into a sep- 
arable equation by changing both the dependent and independent variables if, and only if, 
there exist a nonconstant analytic function W of the complex variable z = x + iy and 
functions $ and ^ such that the canonical invariants and /e satisfy the following 
conditions : 


~ 0 


(308) 


</ E = [$(u) + #(v)] 




dW 


dz 


where 


(309) 


u 


v =j?Ww J 

(C22) If the canonical invariants of equation (31) do satisfy conditions (308) and (309), 
then this equation can always be transformed into a spearable equation by introducing both 
a new dependent variable V defined by 
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V = e w / 2 U 


(310) 


where <o is determined to within an unimportant constant by 


a = u> 


b - co 

y 

and new independent variables | and 77 defined by 


£ = <P( 

?7 = $(v) 




(311) 


where <p and are any convenient nonconstant functions and u and v are determined 
from condition (309) . 

Notice that, by choosing <p(u) = u and $/(v) = v, any equation which satisfies condi- 
tions (308) and (309) can always be transformed into a separable equation by the conformal 
transformation 


z - W 


More generally, the change of variable (311) represents an orthogonal tranformation 
and it is shown in reference 1 that the particular choice of the functions cp and \ p mere- 
ly serves to modify the scale factors of the coordinate geometry associated with the new 
independent variables. 


Direct Calculational Procedure for Testing / E 

As was done in the previous cases, we shall now give an alternate form of condition 
(309) which can be used to test the canonical invariant directly. In order to simplify 
the following analysis , we shall assume that is analytic both in its dependence on x 
and in its dependence on y. Actually, this restriction can be weakened considerably, 
but it is not felt the additional effort is justified. 

Suppose first that satisfies condition (309). Then, since W is nonconstant, 
equation (309) can be written as 
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(312) 


/e 

dW 


2 


= $(u) + 'k(v) 


dz 


But there exist functions $ and ^ 


such that equation (312) holds if, and only if, 



9u 3v 



= 0 


(313) 


We shall be able to greatly simplify the analysis and use many of the results obtained 
in the hyperbolic case if we perform an analytic continuation of the function /e to com- 
plex values of the variables x and y and introduce the complex substitution 


z s x + iy 


z* = x - iy 


(314) 


It is important to notice that z* is the complex conjugate of z only when the vari- 
ables x and y are real. Now when x and y are each allowed to range over the com- 
plex plane, the variables z and z* can certainly be varied independently and we can 
therefore treat them as independent variables. 

Now suppose x and y are real variables and X = h + ig is an analytic function of 
the complex variable z = x + iy. If x* = h - ig is the complex conjugate of X, then 
when the substitution (314) with x and y real is used to eliminate x and y in the par- 
ticular formulas for x and x*, it will be found that the expression for X will contain 
only z and the expression for x* will contain only z*. 

Hence , if the variables x and y are continued analytically into the complex plane , 
then X will be a function of the independent variable z only and X* will be a function 
of the complex variable z* only. However, x* will be the complex conjugate of X only 
when x and y are real, or equivalently, only when z* is the complex conjugate of z. 
If P is any function of the complex variable z, we shall always denote by P* the func- 
tion of z* which is equal to the complex conjugate of P when x and y are real (or 
equivalently, when z* is the complex conjugate of z). 

In view of these remarks, we see that when x and y are extended to complex 
values, then W, dW/dz, etc. are functions of the independent variable z only; and W*, 
(dW/ dz)* = (dW*/dz*), etc. are functions of the independent variable z* only; and, for 
example, (dW*/ dz *) is the complex conjugate of (dW/dz) only when x and y are 
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restricted to the real line. In addition, the relations 


and 


u =lw +iw* 


(315) 


W - W* 
2i 


still hold when x and y are complex. Since W is a non constant function of z only 
and W* is therefore a nonconstant function of z* only, we can also take W and W* 
as independent variables; and equation (315) shows that u and v can be treated as in- 
dependent variables. We shall use the notation W' for dW/ dz and (W*)’ for 
dW */ dz * = (dW/dz)*, etc. 

Now the principle of analytic continuation shows that equation (313) must still hold 
when x and y are continued to complex values. Equation (315) shows that 


_a_ = _a_ t a 
du aw aw* 

dv Law aw*J 


But 


a _ l a 
aw w’ az 


and 


Hence, 


a _ l a 
aw* (w*)» az* 


a 2 __i a_ j_ i a_ a a 

au dv w r az w az (w*)* az* (w*)’ az* 
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r 


Using this in equation (313) we find 


1 9 

J_ _9_ 

/e 

1-1 3 1 

fl 3 

/ E 1 

W' dz i 

|W» 9z 

W'(w*)’ 

j (W*)» 9z* | 

^(W*)’ 9z 

W»(W*)’ 

-*y 


or, since (W*)' is indpendent of z and W’ is independent of z*, 


_9 1_ _9_/j_ 4 \ _ d 1 

dz W* dz \W» E ) dz* (W *)» 

Upon introducing the function S of z only by 


_i L_ / 

dz* (w*)' E 



and according to our convention 



this equation becomes 

23 WL + K W Z + s "s* - 23 *W Z , Z , + 3 <s * ), W 2 . + (s *>” */■ 


(316) 


(317) 


As in the derivation of equation (135) from condition (131), we can show that the 
steps of this argument can be reversed to establish the following conclusions: 

(C23) There exist a nonconstant and analytic function W = u + iv of the complex var- 
iable z = x + iy and functions $ and ’J> such that satisfies condition (309) if, and 
only if, if, there exist a nonzero function S of t he complex variable z only and a non- 
zero function S* of z* only which is equal to the complex conjugate of S whenever z* 
is equal to the complex conjugate of z such that / E satisfies equation (317). 

(C24) If nonz e ro functions S and S* can be found such that equation (317) holds, 
then the analytic functions W of the com plex variable z for which condition (309) is 
satisfied can be calculated from equation (316). 

Thus, if the function S is known, conclusions (C21) and (C22) give a procedure for 
calculating a change in the independent variables which will transform equation (31) into 
a separable equation. 
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Now equation (317) is essentially the same as equation (135) for the hyperbolic case. 
Therefore, all the solutions to equation (317) can be found by the same procedure as was 
used for equation (135). The remarks following conclusion (C9) apply to this case also, 
except that, after the constants have been determined by substituting the solutions into 
equation (317), it must be possible to adjust the constants so that S is not zero and S* 
reduces to the complex conjugate of S when x and y are real. 

Again, the two exceptional forms of the function for which the general proced- 
dure will not work must be considered separately. Thus, when = 0, condition (309) 
can always be satisfied. If 0, we define the function by 


/ e -4- — 4- — / E 

E Se dz I/e 9z * E 


for / E * 0 


or introducing the variables x and y this becomes 



d_ I 1 9 /e\ 
_ 3x Ve Sx / 




(318) 


(319) 


In particular, this shows that like itself the function / E is real whenever the 
variables x and y are real. 

Now suppose there is a constant c Q such that 



(320) 


Clearly, c Q must be real. (Otherwise, the equation could not hold with x and y 
real.) For the present purpose, we need consider only real values of x and y. If 
c Q = 0, then 


Jif JL + ±(-L = o 

3 x\/ e Sx J iy\/ E oy ) 

but this implies that there exists a nonzero analytic function X of the complex variable 
z such that 


/ E = tti1x( z )1 2 
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where 7Tj is a constant. It is clear that condition (309) will always be satisfied if we 
take 


dW 

dz 


X 


If c Q # 0, equations (319) and (320) show that 


9 / 1 

9x \/ E 3x / 




" 4c o<^ E 


It is shown in appendix I that the most general (real) solution to this equation is 

A c o ~ 0 

(321) 

^ / E c o S 0 

where h can be any harmonic function. 

It is easy to see from equation (294) that if we put 


1 < h x > 2+ (\> 2 


2 c 


cosh 2 h 


A 


< 


1 < h x) 2+ < h v) 2 


2c 


cos^ h 


2c q cosh u 


4>(u) = 1 


2c cos u 
V o 


for / E c o - 0 
for <4 c o >0 


and h = u, then equation (321) satisfies condition (309). 

By changing the harmonic function h it is possible to express the solution (321) in 
many different ways. Several of these are given in appendix I. Thus, the second line of 
equation (321) can be written either as 


1 

2c u 2 
o 


dW 


2 


dz 
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or as 


-±[p(u;g 2 ,g 3 ) +p(v;g 2f -g 3 j\ 


dW 

dz 


2 


where W = u + iv is an analytic function of z. It is easy to see that the first of these 

9 

satisfies condition (309) with 4>(u) = l/2c Q u and 'J'(v) = 0, and that the second satisfies 
condition (309) with ^(u) = p(u;g 2 ,g^ )/2c Q and ^(v) = ^0(v;g2, -gg)/2c o . It is shown 
in reference 2 that all the changes of the independent variable which will transform equa- 
tion (31) into a separable equation when 


E 


2c o u 


dW 

dz 


2 


are given by equation (122) which transforms this expression for f E into the form (121). 

We have now shown that 

(C25) The condition (309) is always satisfied if either the invariant or the func- 
tion j E defined by equation ^318) is equal to a constant. 

Now as in the hyperbolic case we again suppose that / E # constant, / E i= constant, 
and that the functions S of z only and S* of z* only are any two simultaneous solu- 
tions of equation (317). The same procedure as was used in the hyperbolic case now leads 
to the equation 


2-^lS + 5kS’ = 2 — S* + 5k*(S*)’ (322) 

dz dz* 


where 






k* / E 

E az* E j 


(323) 


Notice that since / E is real for x and y real, k* becomes the complex conjugate of 
k when x and y are real. This shows that the notation used herein is justified. 

Now if k were zero, equation (323) would show that 
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Hence, we could conclude that for x and y real 

9 /e _ 8 ^E 

3x 3y 


= 0 


or y'g = constant. But this is contrary to hypothesis. Similarly, it can be shown that 
k* =£ 0. Therefore, the procedure used to derive equation (148) can be applied here to 


show that if we define Sq , s j , and 

s 2 by 


( k *)3/5._a_ 

( k *) 2 / 5 / i 


dz* 


s 0 = < 


- (k*) 


&*■>. 

2/5 j 


Jk 


for K = 0 


for K * 0 


3 /5_3_l/k*V 

dz* |\k / K 


k 2 / 5 ! 

L\k*/ J r 


- k for K * 0 


s l = < 


5 1 


2 .2/5 dz 


jL'l k 2/5 
k* 


for K = 0 


(324) 


(325) 


— (k*)*^'* — 

2 dz* 


3 2 " i 


(k *) 2/5 ^ k^ 


K 


for K * 0 


- — for K = 0 
2 k* 


(326) 


with 


K = 



z 


then S must satisfy the equation 
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SqS + SjS' + s 2 S" + 0 


Similarly, S* satisfies the equation 

sgS* + sJ(S*)» + S$(S*)" = 0 (328) 

where the coefficients sg, s|, and s| can be obtained from equations (324) to (326) by 
interchanging the starred and unstarred quantities. 

The solutions to equations (327) and (328) can now be obtained by exactly the proce- 
dure as used for equations (148) and (152). We shall therefore not repeat it here. Since 
an elliptic equation cannot be transformed into an equation which is weakly separable but 
not separable, this completes our discussion of the elliptic case. 


SUMMARY AND CONCLUDING REMARKS 

Necessary and sufficient conditions which a linear second-order partial differential 
equation in two independent variables must satisfy if it can be transformed into a separ- 
able equation (or into an equation which is weakly separable but not separable) have been 
obtained. These conditions together with the appropriate changes of variable (which will 
bring the equation into separable form) are summarized in tables I and II. 

In addition, a procedure has been developed for testing by direct calculation whether 
a given equation of this type can be transformed into a separable equation, and a pro- 
cedure has been given to calculate the new variables. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 14, 1970, 

129-01. 



TABLE I. - TRANSFORMATION OF EQUATIONS INTO SEPARABLE EQUATIONS 


CD 

00 


Type of equation 


Hyperbolic 


Parabolic 


Elliptic 


Canonical form 


3x 2 ay 2 ax ay 


3 u + a— +b — + cU = 0 


3x‘ 


.2 dx dy 


ru + £U 3U 3U 
ax 2 ay 2 


+ a — + b — cU = 0 
3x 3y 


Canonical form 


First 

invariant 


• / H = a y +b x 


y p = b 


• / E= a y- b x 


Second 

invariant 


/ H = C 4 (a x t V‘| (a2 ' b2) 




\ 2b / 



/ E = c -i(a x+ b )-I ( a 2 + b 2) 


Functional forms 
of the invariants 
for equations 
which can be 
transformed 
into separable 
equations a 


First 

invariant 


y H = 0 


P djM x 


Second 

invariant 


/ H = [*<u) + *(v)](«4 - U ^ 


[-L R(<?>)] - 

b / ^y \ 

+ fc 

Lv(y) J x 

2 W_ 

X 


where 


u = 1 F(o) + 1 G(t) 

2 2 

v = 1 F(<j) ~ 1 G(t) 

2 2 

cr = x+y t = x - y 
FW 0 G’ * 0 


where 


dj * 0 


>E = ° 


/ E = [$(u) + ^vjj 


where 

u = i?eW 
v = JtaW 

and where W is any analytic 
function of the complex vari- 
able z = x + iy with 
dW/dZ * 0 


Change of vari- 
able which 
transforms 
differential 
equation into 
a separable 
equation 


Dependent 

variable 


where 


r = e w / 2 U 


x 

b = 


where 


(l/2)[w + (l/2)ln|b| +e] 

V = e U 


1 = f a dx 


r =« -J—R((p) + -fo 

v(y) 2 \(p 


where 


V = e"^ 2 U 


a = w x 
b = w 


, = 1 


ft h*J2c-± 

2 I 2 


a 2 - I — I 



Independent 
variable s a 


£ = $?(u) where £> f * 0 

7] = $(v) where $/’ * 0 


I = 

?7 = 3/(y) where * 0 


£ - where ^ 0 

T] = $(v) where * 0 


a ^, F, G, R, v, dp and (p can be any functions of their arguments. 



TABLE n. - TRANSFORMATION OF EQUATIONS INTO EQUATIONS WHICH 


ARE WEAKLY SEPARABLE BUT NOT SEPARABLE 


Canonical form 


Type of equation, hyperbolic 


i!s. i!s + a iS + b iS + c u = o 


3x 2 3y 2 


9x dy 


Canonical First invariant 
invariants 

Second invariant 


Functional forms of the invariants 
for equation which can be trans- 
formed into an equation which is 
weakly separable but not 
separable 21 


where 


•'h “ a y + b x 


/ H _ c - I ( a x + by) - 1 (a 2 - b' 

y H =$(v)F’(ff)G'(T) 

/ H = *(v)FW(t) 


V = F(cr) + G(t), V = F(ct) , or v = G(r) 


O' = x + y, r = x - y 


F 1 * 0, GW 0, and $ * 0 


Change of variable Dependent 
which transforms variable 
differential equa- 
tion which is 
weakly separable 
but not separable 


V=e IT when v = F + G or when v = F 


J X) = i - i. G'(t) f $(v)dv 
x > d J 


t-0 = 1 r,' 


— + i G'(r) f $(v)dv 

9 A J 


for v = F + G or v = F 


V = e U when v = F + G or when v = G 


F f (cr) J $(v) dv 

f X F'(ff) f <f>(v) dv 
y •). 4 J 


for v = F + G or v = G 


Independent 

variables 3. 


£ = <p(u) where cp' ± 0 

Tj = xj/(y) where 0 


(p ij/ 3>, F, and G can be any functions of their arguments. 


I 


APPENDIX A 

DERIVATION OF EQUATION (140) FROM EQUATION (135) 

For brevity, I will be used in place of /h in this appendix. With this notation 
equation (135) is, after dividing through by I, 



Differentiating both sides of this expression with respect to o and t, we find 



Upon defining p by 



and 
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equation (A2) becomes 


2^££j S + 5 M(j + 4^yj M S’ + 3qS" = 2 ^11 j T + 5 ju t + ju T' + 3juT’ 


Now multiplying equation (Al) by 3jU and subtracting the result from equation (A3), 


we obtain 


) -'Pr '** 


u I 


[j “V/or V 1 


^ = 2„ l(M - 3 m T + sjn - — )t' (A4) 


M I 


Now 
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Interchanging a and r shows that 






= 1 -ifl 4 

M i 4 . 




When these results are substituted into equation (A4) we obtain, after multiplying through 
by I 3 , 



Since we have put I = and p = (I a /I) , definition (136) shows that 



and therefore equation (A5) is the same as equation (140). 
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APPENDIX B 


CONDITIONS FOR WHICH COEFFICIENTS IN EQUATIONS (148) AND (152) VANISH 

Since kj 4- 0 and k 2 * 0, it follows from equation (147) that the condition t^ = 0 
implies that 



* 0 


and that there exists a nonzero function X^ of r only such that 

5 

— = *i( T ) (Bl) 

} 

Hence, it follows from equation (146) and equation (Bl) that the condition t^ = 0 implies 
that 




_d l_/^l\/^2 k 2/^\ 

da ki^ 5 \ k 2/\ k l 2 ) T 


-k 2 = ° 


but this shows that 


_a_ 

da 


L k 2 



- 0 


It follows from this that there exist a nonzero function X 2 of t only and a nonzero func- 
tion y 2 of a only such that 


— = * 2 (T)y 2 (cr) (B2) 

k l 
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and that 



(B3) 


It now follows from equations (145), (B2), and (Bl) that the condition t'^ = 0 implies 
that 




= 0 


Hence, upon using equation (B3), this shows that 


1 


X 





= 0 


This now shows that there exists a nonzero function y ^ of 


o' only such that 


X 

k 


5 

2 

2 

2 




= r 4 ( CT ) 


Hence, upon substituting for kj and kg from equations (Bl) and (B2), we find that 

[M t )| 8 y 4 i«) 

We conclude that there exists a nonzero constant such that 
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and 


x l= x 2»l 


Hence, equation (Bl) becomes 


or using equation (B2) 




5 

ff l X 2 


This equation can also be written as 


_a_ 

9t 









(157) 


But this is essentially Liouville’s equation, and therefore its most general solution is 
(see eq„ (138)) 



\2(r) 




Ag( T )yg(ff) 

Jy 3 (a) - X 3 (r)] : 


(155) 


and equation (B2) now shows that 


k 



1 


X|(T) 

r\(°) 


C X 3 (T)r 3< C,i T 
{ [yjW - x 3< t )] 2 J 


(156) 


Thus, = 0 only if k, and kg satisfy conditions (155) and (156). It is an 

easy matter to verify that the steps of the preceding argument can be reversed to show 
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that equations (155) and (156) are also sufficient conditions. It now follows from the 
symmetry between equations (145) to (148) and equations (149) to (151) and between equa- 
tions (155) and (156) that t^ = t^ = t^ 2 ^ = 0 if and only if, = 0. It is 

also clear that equations (B2) and (157) are necessary and sufficient conditions for the 
coefficients in equations (148) and (152) to vanish. 


101 



APPENDIX C 


CALCULATION OF S AND T WHEN COEFFICIENTS 
OF EQUATIONS (148) AND (150) ALL VANISH 

Upon substituting equations (155) and (156) into equation (142), we get 

r. . 


da |\y, 


X 3 y 3 


L' 


(y 3 - \ 3 y 


dr 


)(X 9 T) 


»8 j l 


(Vq “ X«) 


Carrying out the indicated differentiations and rearranging gives 


y 3' x 3 2 (S_ nl _ 


y 3 


y 3 \ y 2 




— — (TX 2 X' 3 ) + 2(TX 2 X’ 3 ) 


(Cl) 


Differentiating equation (Cl) with respect to both a and r gives 


-X', 


W5V 

_ y 3 \ y 2 >\ 


= y 3 


~ (tx 2 x* 3 ) 


or 


1 


II 

1 

■h 

— (TXoXk)’ 

y 3 

_ y 3 \ y 2 /_ 

X 3 

L 3 , J 


since the left side is a function of a only and the right side is a function of r only there 
exists a constant, say tt 2 , such that 


y 3 


1 W 

L y s U yj 


= Vc 


X 3 


4 (X 2 X 3 T) 
X 3 


TT 


= ~*2 
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Upon integrating these two equations, we get 


Syk 7T n o 

— =Y r 3 +7r 3 y 3 +7r 4 
(™ 2 *3> = X 3 + n 5 X 3 + n 6 


(161) 


where through n g are constants. 

Upon substituting these results back into equation (Cl), we find that 

71 5 = ~ n 3 
= -ttq 


so 


TX 2 X 3 = 


^.2 . 

'Y X 3 " n 3 X 3 


(162) 
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APPENDIX D 


DERIVATION OF EXPRESSION FOR / p 

We first differentiate equation (231) with respect to x to obtain 



Next squaring both sides of equation (231) gives 



Finally, multiplying equation (231) by -b x /2b and using equation (230) gives 



Upon adding these three equations, we obtain 
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Using equation (230) to eliminate cp /cp yields 

Xa X 







+ 






(Dl) 


After differentiating equation (226) with respect to y, we obtain 



II 

v 


<P XV 

+ 2-5Z 



y 


2 

Upon using this equation to eliminate <p/cp„ and equation (226) to eliminate (<p ) in 

xy x x 

equation (Dl), we obtain, after rearranging the left side and dividing by b, 


r * 



Subtracting this result from equation (229) shows that 
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where we have defined the function R of <p only by 



After differentiating this result with respect to x, we obtain, upon recalling that e g 
depends (implicitly through 4>) only on y, 



(D2) 

Differentiating equation (231) with respect to y shows (since dj and d 2 are func- 
tions of <p only) that 



Upon recalling definition (46) and noting that 
tracting equation (D3) from equation (D2), that 


by/b) = (k> x /b) y > we find > afte 1, sub_ 



which is the desired result. 
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APPENDIX E 


DERIVATION OF ALTERNATE FORM OF CONDITION (233) 

Suppose first that <p v * 0. Now there exists a function of cp only such that 
equation (233) holds if and only if , 


9 / b 


9y l 2 

V^x 


a / b 


<P.~ ^ =° 




Upon differentiating by parts, shows that this is, in turn, equivalent to 


- 0 

dy\w x ) dx \hcpl 7 


Another differentiation by parts, we find that this is equivalent to 


i. b (h. 


ax \<p 


^v\ a / b 


<p x ) 9 y \ v 


Hence, we can conclude from this that, if cp^ + 0, there exists a function d^ of cp 
only such that equation (233) holds if and only if, 


b (h\ .li JS\.» T (s 


2 \v 


Now suppose that cp = 0. This implies that cp is a. function of x only. If there exists 
a function d^ such that equation (233) holds, then this implies (since d^ depends on x 
or y only implicitly through cp) that d/v is a function of x only. Hence, 


But this implies that equation (El) holds. Conversely, if equation (El) holds, then the 
fact that <py = 0 shows that 


■ ■■ ■ II ■ ■■ ■ ■■■■■■■■■■ ■ m mi ■ 



b 

v 


= 0 


since v + 0„ Thus, b/v and cp are functions of x only. Hence, 


b 



is a function of x only and therefore (since <p is a function of x only) a function of (p 
only. This shows that there exists a function dj of cp only such that equation (233) 
holds. We can therefore conclude that there exists a function d^ of cp only such that 
equation (233) holds if and only if, equation (El) is satisfied. 
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APPENDIX F 


APPLICATION OF TRANSFORMATION (240) TO EQUATION (30) 

It follows from equations (41) to (43) and (234) that under the change of variable (240) 
the coefficients a, la, and c of equation (40) are given by 


~ / 1 b x^ 

a = - w Y + 

V x 2b; 


(FI) 


E = b 


(F2) 


and 


2c = 2c - a f u) + — — —■ + a\ 
\ 2 b / 


- b Kif +9 y 


T +a ) x + K“ x+ ^ ta ) 


= 2c +1 (u> +1 — - - a 2 - b fu > +1 -1 + e V (w +- — + a^ 

2 \ X 2b/ 2 V y 2 b y /V X 


2 b 


x 


2c - - 


a 2 - 

\2b , 


L +1 ^ - b (e„ + +— +- — o>„ - bG0 +1 b„ - w. 


2 b 


x 


2 2 b 


y xx 


Upon substituting in definition (235) this becomes 


■ b 1 

2 c = b ( 52 + — - co \ - I co - — — — co b 

y) v XX 2 b x 5 y l 


2b 


Hence, equations (238) and (239) now show that 


~ b 
c = — 


2v 


R (cp) + — v’ 
2 


(F3) 


Substituting equations (Fl) to (F3) into equation (40) yields 
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v YV - U x +- — Vv +bv v + — [*%) +- v’i v = o 
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APPENDIX G 


DERIVATION OF EQUATION (251) FROM EQUATION (249) 

Upon differentiating equation (249) with respect to x, we get 

( V ^ p ) xy " bT yyx " b x T yy ' b y T xy ' b xy T y = (Vp)^ 


This equation can be rearranged to obtain 



Hence, 



Substituting equation (250) gives 
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or 



Differentiating equation (250) with respect to x yields 



Substituting this result into equation (Gl) gives 
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Now 



Hence, if we define / p by 
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l 




equation (G3) becomes 


2v’/ p 



(252) 


(251) 
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APPENDIX H 


PROOF THAT EQUATIONS (249) AND (250) CAN ALWAYS BE SATISFIED IF(/ p ) = 0 

x 


Let the domain of definition D of equations (249) and (250) be divided into two sub- 
domains D + and D” such that 

b(x,y) > 0 for (x,y) e D + 


and 


b(x,y) < 0 for (x,y) e D" 

We shall show that if = 0, equations (249) and (250) always possess nonzero 

solutions which contain three arbitrary constants for (x,y)eD + . A similar proof will 
show that they also possess nonzero solutions for (x,y) e D”. We can conclude from this 
that these equations possess nonzero solutions in D which involve three arbitrary con- 
stants but A\hich may have to be specified by different functions in different parts of D. 
Hence, suppose that (x,y)^D + and (j-pj = 0. Then there exists a function y of 

y only such that /p = y. 

Both equations (251) and (257) now reduce to the same equation, namely, 


2v'y + vy' +— v ,M = 0 
2 


(HI) 


Now the derivation of equation (251) from equations (249) and (250) shows that for 
every solution T of equation (250) (and this equation certainly always possesses a solu- 
tion when v is any given function of y) 


_d_ 

i3x 


1 ^ 

2 b) 


(r/p). - K)„ - (Vp), 


= Vb 


ax 


-h S V/p) y ' G 


bT y) y - ( T >p), 


= 2v'y + vy' +— v’** 
2 


(H2) 


since b ^ 0. 

Now equation (HI) shows that equation (257) has infinitely many nonzero solutions. 
(Recall that a third-order linear ordinary differential equation has three linearly inde- 


115 



Ill I 


Fiji 


pendent solutions and that any linear combination of these solutions is also a solution.) 
Hence, let the function v Q of y only be any solution to equation (257). 

Then equation (H2) shows that there exists a function r T of y only such that 


(Vp) v - ( bT y) v - (Vp) x = r T<y> V5 


(H3) 


for each solution T of the equation 


T +l^T=lv i-iv' 

x bb 2 ° b 2° 


(H4) 


(This equation certainly has a solution T.) 
It follows from definition (252) that 


r(y) = 


Vb 


*7= Sp} --Hf 
vVb P / x 2VVby 


yj 


(H5) 


Now let the function p of y only be any solution of the ordinary differential equation 


Pyy +rP = r T 


(H6) 


and define the function T^ by 


<p(0) = >p + p(y) 

Vb 


(H7) 


Hence, the function T'®' can involve two arbitrary constants. 
Substituting equation (H7) into equation (H4) shows that 


T (0) + 1 ^x T (0) = 1 v ^y _ 1 v , 


2 b 


2 ° b 2 ° 


(H8) 


Substituting equation (H7) into equation (H3) shows that 


( V o*l - 


bT 


(0) 


T (0) / P | 


■Jx 


= r T (y) Vb - 


b -L 

Vv^J 




(H9) 
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Using equation (H5) to eliminate 


(/p/ Vb) x 


in equation (H9) gives 



Vb ( r T 


pr) - 



1 

2 



P 


- V5 (r T - pri - (* Py -iip) y -ip(^ 

= Vb (r T - py - p yy ) 


Hence, equation (H6) now shows that 





(H10) 


Thus, we have shown that, if y/'-pj = then for each solution v Q of the ordinary 

differential equation (257) there exists a function (containing two arbitrary constants) 

such that v Q and satisfy equations (249) and (250). Notice that this proof was con- 
structive and can therefore be used in practice to find the function once a solution 

v^ to differential equation (257) is known. 
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APPENDIX I 


SOLUTIONS TO ELLIPTIC LIOUVILLE'S EQUATION 

The elliptic form of Liouville's equation is 

±(± + ±(1 9U\ = c TJ 

3x \U ax/ 3yVU dyj 


(II) 


where c Q is a constant. Let 


X 1 =h 1 + ig : 


be any nonconstant analytic function of the complex variable 

z = x + iy 

and let H be a nonzero function of hj and gj. Put 


U = H 



( 12 ) 


Then since In | dXj/dz | is a harmonic function, 

~_d_ (1 J!*L\ - (1 8H \ 

ah^H 3hJ + a gl \H a g J 



Hence, equation (II) becomes 


JL/I 5H \ + _L (1 

dh^E 9hJ dg 1 \h a g J 


C 0 H 


(13) 


If we put H = F, where F is a function of hj only, then equation (13) becomes 


F" 
. F 



F 


2 


c o F 


(14) 
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Hence, upon setting F' = P, we find 


F M = P — 
dF 


and equation (14) becomes 


1 1 d(P 2 ) P 2 _ c F 

2 F dF f 2 ° 


or 


(F*) 2 = P 2 = F 2 (2c d F + Cj) 


where c ^ is a constant of integration. 
This shows that 


± 



+ c 


2 


Thus, 


Hence , 



if Cj + 0 



2 (ftf 

1 

F 

C ° ' 2 ' cosh 2 

< h l + c 2> 


if Cj * 0 


2 

. c o (h l + c 2 )2 


if Cj = 0 
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The solution (12) of equation (II) is therefore 


r 


2 _ 

c„ 





V 2 ) 

2 

dXj 

dz 

2 

cosh^ 

< h l + c 2>| 


if Cj + 0 


dX„ 


dz 


if = 0 


C 0 (hj + c 2 r 

Now put Cg = | Cj | and define the harmonic conjugate functions h and g by 


h = 


(hj^ + c 2 ) if Cj * 0 


hj + c 2 if Cj = 0 


r 


g = < 


2 


gj if Cj * 0 


g 1 if c ! = 0 


If we define the analytic function X of the complex variable z by 

X = h + ig 


the solution (15) becomes 


2 I dz | 
c o cosh 2 h 


if Cj > 0 


if Cj = 0 


2 I dz I 
c o cos 2 h 


if c x < 0 


for any nonconstant analytic function X. Actually, since the analytic function X is quite 
arbitrary, the expressions on the second and third lines of (16) are the same, for if we 
put X = 1/i lnXg, hg = Q.&X 2 where X£ is any analytic function. Then, 


X 2 = e l 


2 = . e ix dX 
dz dz 


dX 2 = e i(x-X*) IdX I 2 


cos 2 h = = (e 


i[(X+X*)/2] _-i[(X+X*)/2]'l 2 


i i(X*-X) /e iX + e" ix * 


_ e i(X*-X) / X 2 +X 2 


= e^ x * -x ^(h 9 ) 2 
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Hence, 


dx 2 

2 

dX 

dz 

_ 

dz 


(h 2 ) 2 cos 2 h 


Notice that when the solution U is given by the first line of equation (16) 

UCqSO 

and when the solution U is given by the second (or third) line of equation (16) 

Uc 0 ~ 0 


Thus, these forms of the solution cannot be transformed into one another simply by 
changing the arbitrary analytic function X, and they therefore represent totally different 
expressions for the solution. Hence, 


r 


dX 

dz 


c o cosh 2 h 


U = < 



dX 

2 

2 

dz 



L o cos h 


if Uc Q ^ 0 


(17) 


if Uc Q ^ 0 


for any analytic function X. 

We shall give an alternate derivation of equation (17) and at the same time demon- 
strate that this is the most general real solution to equation (II). It is shown in refer- 
ence 3 (p. 194) that the general solution to equation (II) is 

dX dA 

U = — — dz * (I 

C o (X + A) 2 

where X = h + g is any nonconstant analytic function of the complex variable z and A 
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is any analytic function of the complex conjugate variable z* = x - iy. This solution is 
not real valued, however, unless the functions X and A are related. In order to find 
this relation, notice that for any function A there exists an analytic function T of the 
complex conjugate variable X* such that 

r[x*(z*)] = r(z*) 


Hence, equation (18) becomes 


dx dX* r , 

dX 

T' 

8 dz dz* 

dz I 



'o (x + ry 


(X + ry 


(19) 


where 


= dr 

dX* 


* 0 


Now U is real if, and only if, 


r» (ry* 

(x + r ) 2 (x* + r *) 2 

or 

x + r _ x* + r* 

( F ’) 1 / 2 ( p * 1 )^/ 2 

Now this equation can be extend to complex values of h and g, and the variables X and 
X* are then independent. Recalling that T* and (T 1 )* are functions of X only, we 
find, upon differentiating equation (I- 11) with respect to X*, that 


( 110 ) 


( 111 ) 


r 

t 

+ X 

1 

f 

1 

^r*) 1 / 2 

LnH 

(r**) 1 / 2 


Upon differentiating again with respect to X* we obtain 
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r 

(Ft ) 1 / 2 

Now if [l/(r’) 1//2 ] M * 0, this equation would show, upon division by [l/(r’)^^ 2 ] , that 
X was equal to a constant. Since X is a nonconstant function, we conclude that 


+ X 


V?v 


= 0 


Ur») 


nl/2 


= 0 


Upon integrating this equation, we find that 


r 


r = < 


CjX* + c 2 


+ Cg if C-^ # 0 


c^X* + Cg if Cj = 0 


( 112 ) 


where Cj to Cg are constants of integration. Upon substituting these results into equa- 
tion (111) we find that 


(c 3 +X)(CjX* + c 2 ) + 1 



(c| + X*)(cpi + c|) + 1 



(113) 


and 


X + c 4 X* + c 5 X* + c|X + c 



Differentiating equation (113) with respect to X shown that 

CjX* + c 2 c*(X* + c^) 

1 Wi 


(114) 


Hence, Cj is real and c 2 = Cjc|. 

Equation (114) shows that there exists a real constant Cg such that 
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C 4 


-i2c 


= e 


6 


and a real constant such that 


-ic 


c 


5 " e 



Hence, equation (112) becomes 


c 3 + 


r = < 


CjCX + c 3 ) 


- if c, ^ 0 
* 1 


“i2ce -2ic fi 

e a* + 2c^e if Cj = 0 


where c^, Cg, and are real. 

Upon substituting equation (115) into equation (19), we obtain 


8c* 


dX 

dz 


if c, # 0 




' 0 (cJX + Cgl 2 +l) 

,2 


dX 


dz 


ic fi \ / ic fi 

e X + CrjJ + ( e X + Crj 


if c^ = 0 


Hence, upon defining the analytic function Xj = h^ + ig^ by 


n 


c* (X + Co) if c* * 0 


Xi = < 


ic fi 

e D X + c^ if Cj = 0 


(115) 


we obtain 
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mill 


Ffl 


IdX, 


dz 


-° (iXjl 2 ^!)' 

IdX, 


if Cj > 0 


U = < 


dz 


if Cj < 0 


0 ( l ■ l x il 2 ) 


dX, 


dz 


“o (hj)' 


if Cj = 0 


(116) 


Notice that, for the first form of the solution (116), Uc Q ^ 0; and for the second two 
forms, Uc Q ^ 0. In fact, these latter two forms are equivalent since they can be trans- 
formed into one another by the conformal transformation 


x 2 - 


1 - X 

1 +X 


1 

2 


Hence, 


r 


8 







C o (hj) 2 


if Uc Q < 0 


if Uc s o 
O/ 


(117) 


We have already shown that the second form of the solution (117) is the same as the 
second form of the solution (17). We shall now show that the first form of the solution 
(117) (which is the form given in ref. 13) is the same as the first form of the solution (17). 
To this end put 


126 



where X = h + ig, and obtain 


X = In X 1 


4 4 4 e -(X.+X*) 

(iXjl^l) 2 (e x e X %l f e X+X* j e (X+X*)/2 + e -[(X+X*)/ 2 ]J^ cosh 2 h 


dx i| _ x+x* 

dX 

dz 1 

dz 


which proves the assertion. 

Hence, we conclude that the most general real solution to equation (II) is given by 
equation (17), with X = h + ig an arbitrary analytic function of z. 

There are many different expressions which can be given for the solution (17). These 
can be obtained by varying the arbitrary function X. For our purposes, the most useful 
of these, which is an alternate for the second form of equation (17), is the one which is 
analogous to that obtained in the hyperbolic case. To obtain this expression, it is 
more convenient t6 start with the forin of the second part of equation (17) given in the cen- 
ter of equation (16), which is 


or, equivalently. 




I dz I 

(X + X*) 2 


(118) 


(119) 


We now define the analytic function W = u + iv of the complex variable z by 



( 120 ) 
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where g 2 and g 3 are real constants and the path of integration maj 
does not pass through a zero of 4t 3 - g 2 t - g^. Then (ref. 12) 

ix = ^(i w ;g 2 ,g 3 ) 

Now it is shown in reference 12 that jp( Z) is real when Z is real. 

[p( z)]* = P(z*) 

Hence, 

-IX* =J3(1»»;8 2 ,8j) 

Upon using these results in equation (119), we find 


w* 


U = - — 


dW 


12 


,<Kl w )-<Kl w * 

The results of example 1 on page 456 of reference 12 now show that 

ki w 4 w *Mi w -i 


u = 


w*' 


jl 

dW 

l\ 

dz 


= — [^( u ;g 2 > g 3) - <P( iv ;g2’ g 3 } ] 
c o 

But example 2 on page 439 of reference 12 shows that 


dW 


dz 


p{ iv;g 2 ,g 3 ) = -p(v;g 2 ,-g 3 ) 


Hence , 


be any curve which 


This shows that 
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U = — [^(u;g 2 ,gg) +^(v;g 2 ,-g 3 )] |^| (121) 

c o 

This form of the solution to equation (II) could also have been obtained by choosing the 
function H to be of the form 


H(h 1>gl ) = F(h 1 ) +G( gl ) 

It is not hard to show by using the results on page 453 of reference 12 that, if we had de- 
fined the function W by 


W(z) s 2 


/ +ix(z) 

* 


dt 


apt^ + 4ajt^ + 6a 2 t^ + 4a^t + a^ 


( 122 ) 


(where a Q to a 4 are any real constants) instead of by equation (120), then there would 
still exist real constants g 2 and g 3 such that equation (121) holds. 
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APPENDIX J 


SYMBOLS 

A coefficient of first derivative in general partial differential equation 

a coefficient of first derivative in canonical partial differential equation 

a transformed coefficient of first derivative in canonical partial differential 
equation 

B coefficient of first derivative in general partial differential equation 

b coefficient of first derivative in canonical partial differential equation 

b transformed coefficient of first derivative in canonical partial differential 
equation 

C coefficient of linear term in general partial differential equation 

c coefficient of linear term in canonical partial differential equation 

c transformed coefficient of linear term in canonical partial differential equation 
c Q constant 

D domain of definition of differential equation 

d r functions of £ for r = 1, 2, 3, . . . 

e r functions of 77 for r = 1 , 2 , 3 , . . . 

F function of a 

f function of | and 77 

G function of t 

g JWx 

g r constants, r = 1, 2 

H function of 77 only in trial solution 

h /?eX 

I special notation for canonical invariant used in appendix A 

i V-i 

J canonical invariant 

J canonical invariant 

j index in operator (can take on values -1, 1, 0) 
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quantity defined in terms of canonical invariant 

quantity defined in terms of canonical invariant 

quantity defined in terms of canonical invariant 

quantity defined in terms of canonical invariant for r for r = 1, 2 

linear operator defined by eq. (54) 

constant 

function of u for r = 1, 2 

Weierstrass ^-function 

function of v for r = 1, 2 

function of a, function of complex variable z 

coefficients in ordinary differential equations for r = 0, 1, 2 

function of r 

variable of integration 

coefficients in ordinary differential equations for r = 0, 1, 2 

dependent variable 

transformed independent variable 

transformed dependent variable 

transformed independent variable 

analytic function of the complex variable z 

function of x and y defined in terms of u and v by eq. (193) 
independent variable in canonical partial differential equation 
coefficients in ordinary differential equation for r = 0, 1, 3 
independent variable in canonical partial differential equation 
complex variable, x + iy 

coefficient of second derivative in general partial differential equation 

coefficient of mixed derivative in general partial differential equation 

coefficient of second derivative in general partial differential equation - general 
purpose function of ex 

independent variable in general partial differential equation 



X function for transforming dependent variable - general purpose function of r 

fi special symbol used in appendix A 

| independent variable in general partial differential equation 

E function of ij in trial solution 
n r constants for r = 1, 2, 3 

a x + y 

t x - y 

$ function u 

cp function connecting £ with x and y 

X function of the complex variable z 
function of v 

xp function connecting r\ with x and y 

O defined by J-p = 3f2/ 3x 

quantities depending on invariant f for r = 1, 2 
co function used in transforming dependent variable 

Subscripts: 

E elliptic 

H hyperbolic 

P parabolic 

Superscripts: 

* complex conjugate or quantity which reduces to complex conjugate when independ- 
ent variables x and y are real 

’ denotes differentiation with respect to argument (prime) 
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